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ABSTRACT.  The  dynamic  behavior  of  the  mercury  damper 
is  investigated.  Particular  attention  is  paid  to  the 
eccentric  annular  mercury  configuration,  which  is  the 
final  continuous  ring  phase  which  occurs  in  the  opera¬ 
tion  of  -rll  mercury  dampers.  In  this  phase  the  damp¬ 
ing  is  the  poorest,  and  the  system  is  closely  linear. 

The  crescent,  or  broken-ring,  case  is  also  considered. 

During  the  course  of  the  investigation  the  hydro- 
dynamic  problem  is  treated  as  three  dimensional,  and 
extensive  use  is  made  of  a  variational  principle  of 
least  viscous  frictional  power  loss,  which  is  derived. 

A  variational  principle  of  least  constraint  is  also 
used  to  advantage.  Formulas  for  calculating  the  be¬ 
havior  of  the  mercury  damper  are  obtained. 
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INTRODUCTION 


The  purpose  of  the  present  investigation  is  to  examine 
the  dynamic  behavior  of  the  mercury  damper,  and  to  devise 
means  for  calculating  this  behavior.  Since  mercury  dampers 
can  operate  in  several  different  ways,  depending  upon  the 
mercury  configuration,  the  nature  of  the  inner  free  surface, 
and  the  state  of  turbulence;  and  since  any  given  damper 
operates  in  at  least  two  very  different  manners,  the  over¬ 
all  behavior  of  a  mercury  damper  may  appear  to  be  both 
complex  and  peculiar. 

It  will  be  seen  that  any  mercury  damper  acts  more  ef¬ 
fectively  when  the  mercury  forms  a  broken  ring  than  when  it 
forms  a  continuous  ring.  It  follows  that  the  transition 
from  a  broken  ring  to  a  continuous  ring,  which  occurs  in  any 
case  during  the  course  of  the  damping,  corresponds  to  a 
transition  from  a  more  effective  to  a  less  effective  action 
of  the  damper.  It  is  particularly  important  that  the  be¬ 
havior  of  the  damper  in  this  unfavorable  "continuous  ring" 
case  be  analyzed;  and  a  considerable  portion  of  this  report 
is  devoted  to  this  analysis.  The  more  favorable  broken  ring 
case  is  also  treated;  however,  it  is  very  much  simpler  than 
the  unfavorable  case  just  mentioned. 

It  will  be  shown  that  the  viscous  frictional  power  loss 
in  the  mercury  is  of  primary  importance  in  eliminating  the 
precession,  or  wobble;  also  that  this  loss  satisfies  a  mini¬ 
mum  principle  which  constitutes  a  powerful  tool  for  deter¬ 
mining  the  velocity  distribution  in  the  mercury,  and  hence 
the  power  loss  itself.  The  pressure  distribution  is  deter¬ 
mined  using  another  variational  principle--the  principle  of 
least  constraint--which  is  devised  for  this  purpose. 

It  will  be  shown  that  when  the  wobble  is  small  and  the 
mercury  is  in  the  unfavorable  continuous  ring  configuration 
the  system  is  approximately  linear,  and  its  behavior  is 
described  by  linear  differential  equations  with  constant  co¬ 
efficients.  The  usual  methods  for  handling  such  equations 
can  hence  be  applied.  This  is  not  true  when  the  wobble  is 
larger  and  the  mercury  is  in  the  broken  ring  configuration. 
In  both  cases,  however,  specific  formulas  are  obtained  for 
calculating  the  behavior  of  the  system  (see  Appendix  B). 
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MECHANICAL  SYSTEM.  STATEMENT  OF  PROBLEM 


The  mercury  damper  consists  of  a  circular  channel,  or 
race,  partially  filled  with  mercury,  as  shown  in  the  figure 
below.  Two  such  dampers  are  mounted  with  their  planes 


mercury 


perpendicular  to  the  axis  of  symmetry  of  the  main  spinning 
body,  and  with  their  centers  lying  on  this  axis  on  opposite 
sides  of  the  center  of  gravity  of  this  body  and  equally  spac¬ 
ed  from  it,  as  shown  in  the  following  figure.  By  using  two 


2 


NAVWEPS  REPORT  8611 


dampers  the  translational  effects  on  the  main  spinning  body 
cancel,  whereas  the  rotational  effects  add.  We  shall  suppose 
that  the  ratio  of  the  principal  moments  of  inertia  of  the 
main  body  at  its  center  of  gravity  and  the  initial  motion 
given  to  this  body  are  such  that  this  body  precesses  about 
an  axis  which  passes  through  the  center  of  gravity  and  ap¬ 
proximates  the  spin  axis,  the  angular  velocity  of  precession 
exceeding  that  of  spin.  More  exactly  let 

A  =  moment  of  inertia  of  the  main 

body  about  its  axis  of  symmetry, 

B  ~  moment  of  inertia  of  the  main  body 

about  a  transverse  axis  through  its  (1) 

center  of  gravity, 

v  =  angular  velocity  of  spin, 

P  =  angular  velocity  of  precession; 

then 

-g  >  1  and  p  =  -^g-  closely.  (2) 

It  will  evidently  suffice  to  consider  one  of  the  two  dampers; 
hence  looking  down  upon  the  upper  damper  shown  in  the  figure 
at  the  bottom  of  page  2  the  channel,  or  race,  containing  the 
mercury  appears  as  shown  below.  The  center  of  the  channel 
is  S,  which  point  is  also  the  point  of  intersection  of  the 
spin  axis  with  the  plane  of  the  damper.  The  channel  rotates 


channel 
or  mercury 
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about  point  S  with  angular  velocity  v;  however,  in  addition 
point  S  rotates  about  W,  the  point  of  intersection  of  the 
precession  axis  with  the  plane  of  the  damper,  with  angular 
velocity  /3.  We  have  seen  that  /?  >  v . 

Let  us  observe  the  damper  from  a  reference  system  which 
rotates  about  the  precession  axis  with  angular  velocity  13 
counterclockwise;  then  the  velocity  vector  V  and  the  acceler¬ 
ation  vector  A  of  any  particle  of  mercury  will  be  composed 
of  components  as  follows: 

V  -  VR  +  VL> 

A  =  *R  +  AL  +  AC 

where 

V  and  A_  are  the  relative  velocity  and 
K  relative  acceleration  vectors, 

respectively,  seen  while  riding 
with  the  moving  reference  system, 

Vy  and  A  are  th^  locked  velocity  and 

locked  acceleration  vectors,  re¬ 
spectively,  which  the  particle 
would  have  if  it  were  locked  with 
the  moving  reference  system, 

A^  ■-  2/f  x  VR  =  Coriolis  acceleration  vector 

P  -  precessional  angular  velocity  vector. 

The  acceleration  components  Al  and  A q  can  evidently  be  re¬ 
placed  in  effect  by  two  systems  of  inertia  forces  of  inten¬ 
sities  (-Al)  and  ( - Aq )  per  unit  mass,  respectively.  Our 
problem  thus  becomes  one  in  which  we  may  regard  ourselves  as 
being  stationary,  and  the  velocity  and  acceleration  of  a 
mercury  particle  as  being  given  by  Vr  and  Ar ,  respectively. 
The  main  body  now  apparently  rotates  clockwise  about  the  spin 
axis,  which  passes  through  the  stationary  point  S,  with  an 
angular  velocity  (/J-i/).  Two  body  forces  are  acting:  the 
centrifugal  force  of  intensity  (-Al)  per  unit  mass,  and  a 
"Coriolis  force"  of  intensity  (-A^). 

The  centrifugal  fore*,  which  corresponds  to  the  locked 
acceleration,  is 

Centrifugal  force  *  r/31 1  per  unit  mass  (6) 


(3) 

(4) 


(5) 
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where  r  is  the  radial  distance  from  W,  and  lr  is  a  unit  vector 
pointing  radially  outward  from  W.  Equation  6  is  evidently  the 
negative  gradient  of  the  potential  function 

♦  =  ■  2  r  0  per  unit  mass.  (7) 

The  Coriolis  force,  which  corresponds  to  the  Coriolis  acceler¬ 
ation,  is 

Coriolis  force  =  -  20  x  VR  per  unit  mass.  (8) 

Let  us  temporarily  neglect  the  friction  between  the 
mercury  and  the  channel;  then  the  mercury  will  appear  to  be 
at  rest,  and  VR  =  0,  AR  =  0.  The  Coriolis  force  field  there¬ 
fore  vanishes,  and  the  mercury  is  at  rest  under  the  influence 
of  the  centrifugal  force  field  alone.  Since  the  force  acting 
on  a  particle  of  mercury  which  lies  on  the  free  surface  must 
act  normal  t:o  this  surface  it  follows  from  Eq.  6,  or  from 
Eq.  7  and  the  fact  that  the  surface  must  be  one  of  constant 
<f>,  that  r  is  constant  on  this  surface.  Seen  from  above,  the 
free  surface  thus  appears  to  be  circular.  For  sufficiently 
small  values  of  the  distance  WS-  the  mercury  is  shaped  like 
an  eccentric  cam,  the  outer  circular  boundary  being  centered 
on  S,  and  the  inner  circular  boundary  being  centered  on  W  as 
shown  in  the  figure  below.  The  outer  boundary  is  the  channel 
wall;  the  inner  boundary  is  a  free  surface,  the  radius  of 
which  is  dictated  by  the  amount  of  mercury  in  the  damper. 


channel 

wall 


r 
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As  WS  is  increased  a  situation  arises  wherein  the  two 
circles  become  tangent,  after  which  the  mercury  becomes 
ere  scent- shape'5 ,  Since  the  total  volume  of  mercury  remains 
constant,  rhe  radius  of  the  inner  circle  increases  after  the 
mercury  has  become  crescent-shaped. 

Other  mercury  configurations  will  be  obtained  if  the 
inner  radius  of  the  channel  is  so  large  that  it  forms  part 
of  the  i  •\er  boundary  of  the  mercury  surface,  which  could 
happen  in  the  eccentric  annular  shape,  the  crescent  shape, 
or  both.  We  thus  have  four  possible  cases,  which  arise  from 
the  fact  that  for  either  the  eccentric  annular  or  the  cres¬ 
cent  configuration  the  inner  surface  of  the  mercury  can  be 
either  entirely  free  or  only  partially  free. 

Next  let  us  suppose  that  the  friction  between  the  mercury 
and  the  channel,  which  was  neglected  above,  is  now  "gradually 
turned  on;"  then  since  the  channel  is  apparently  rotating 
clockwise  with  an  angular  velocity  (P-v)  the  mercury  tends  to 
be  dragged  clockwise.  If  the  mercury  is  in  a  crescent  con¬ 
figuration  the  channel  rotates  clockwise  relative  to  the 
mercury  body  as  a  whole  with  an  angular  velocity  (P-v),  and 
the  center  of  gravity  of  the  mercury  will  shift  somewhat  off 
of  the  extension  of  line  WS  in  the  clockwise  sense. 

^f  the  mercury  is  in  an  eccentric  annular  configuration 
t1  e  mercury  flows  in  the  clockwise  sense,  the  result  being 
to  reduce  the  velocity  of  the  mercury  relative  to  the  channel. 
The  mercury  is  dragged  clockwise,  and  its  axis  of  symmetry 
shifts  from  the  line  WS  in  the  clockwise  sense.  Supposing 
that  the  inner  free  surface  of  the  mercury  can  still  be  repre¬ 
sented  with  sufficient  accuracy  by  a  circle,  it  follows  that 
I,  the  center  of  this  circle,  shifts  from  W  to  a  new  position 
which  lies  somewhat  above  the  line  WS. 

The  various  mercury  configurations  described  above  in  the 
case  of  no  friction  are  still  possible  with  friction  acting. 

In  addition  the  mercury  flow  can  in  each  case  be  either  lami¬ 
nar  or  turbulent.  In  view  of  this  and  the  fact  rhat  the  mer¬ 
cury  may  be  either  eccentric  annular  or  crescent  shaped,  and 
the  inner  mercury  surface  may  be  either  entirely  free  or  only 
partially  free,  we  see  that  the  mathematical  analysis  of  the 
mercury  damper  involves  the  consideration  of  eight  distinct 
cases.  Which  of  these  cases  will  exist  at  any  time  depends 
upon  v,  A/B,  WS,  the  inner  and  outer  radii  of  the  channel, 
and  the  volume  of  mercury  in  the  damper.  In  any  given  design 
A/B,  the  channel  dimensions,  and  the  volume  of  mercury  are 
fixed;  also  v  is  approximately  fixed.  As  the  damper  functions, 
however,  WS  decreases  from  its  initial  value  to  zero,  and  the 
velocity  of  the  mercury  relative  to  the  channel  drops  to  zero. 
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During  this  time  all  eight  of  the  above  cases  will  not  be 
experienced,  but  three  may  be.  In  all  damper  designs  the 
last  case  will  be  that  in  which  the  mercury  shape  is  eccen¬ 
tric  annular,  the  inner  mercury  surface  entirely  free,  and 
the  flow  laminar.  Preceding  this  may  be  one  wherein  the 
mercury  shape  is  eccentric  annular,  the  inner  mercury  sur¬ 
face  either  entirely  or  partially  free,  and  the  flow  turbu¬ 
lent.  Finally,  preceding  this  may  be  one  wherein  the  mercury 
is  crescent-shaped,  the  inner  surface  but  partially  free,  and 
the  flow  turbulent. 


ACTION  OF  THE  DAMPER  UPON  THE  MAIN  SPINNING  BODY 


We  shall  now  calculate  the  force  and  torque  exerted  by 
the  mercury  upon  the  channel  walls,  and  hence  upon  the  main 
spinning  body.  In  so  doing  we  shall  view  the  system  while 
rotating  counterclockwise  with  angular  velocity  0,  as  de¬ 
scribed  in  the  previous  section.  We  thus  see  the  relative 
velocity  Vr  and  the  relative  acceleration  Ar;  and  have  a 
centrifugal  force  field  given  by  Eqs.  6  or  7,  and  a  Coriolis 
force  field  given  by  Eq.  8.  Let 

F  =  total  force  exerted  by  the  channel  walls 
upon  the  mercu:  y, 

Fl  *  total  centrifugal  force  exerted  upon  the 
mercury, 

Fc  =  total  Coriolis  force  exerted  upon  the 

mercury,  (9) 

VRq  *  relative  velocity  of  the  center  of 
gravity  of  the  mercury, 

A_,0  =  relative  acceleration  of  the  center  of 
gravity  of  the  mercury, 

M  *  total  mass  of  mercury; 

then  for  translational  equilibrium  of  the  mercury  mass  taken 
as  a  whole  we  have 

F  +  fl  +  Fc  -  M  *w>-  <10> 

Although  we  may  see  the  mercury  as  flowing  in  some  way,  its 
center  of  gravity  will  appear  to  be  fixed;  hence  A^0  «*  0  and 
Eq.  10  becomes 

F  "  “  fl  “  FC*  (11) 
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Multiplying  Eq.  6  by  dm,  the  element  of  mass,  and  inte¬ 
grating  over  the  entire  mass  of  mercury,  we  obtain1 

Fl  =  ?  r  lr  dm.  (12) 

But  rQ,  the  vector  from  W  to  the  center  of  gravity  of  the 
mercury  is  given  by 


"  S  IM  r  3r  d”- 

Equation  13  in  Eq.  12  gives 

Fl  =  M  0*  r0,  (14) 

which  is  identical  with  the  centrifugal  force  that  would  be 
obtained  if  the  entire  mass  were  concentrated  at  the  center 
of  gravity. 

Turning  to  the  Coriolis  force  Fc  we  obtain  from  Eq.  8 


1 


/3  x  VR  dm 


dm. 


(15) 


But,  noting  Eq.  9, 


IMV- 


=  total  momentum  of  mercury  =  M 


=  0, 


(16) 


since  Vrq  *s  0  due  to  the  fact  that  the  center  of  gravity  of 
the  mercury  appears  to  be  stationary.  Equation  16  in  Eq.  15 
gives 


Fc  *  0. 

Equations  14  and  17  in  Eq.  11  now  give  finally 

F  «  -  M 


(17) 


(18) 


1  See  Appendix  A. 
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Next,  we  shall  determine  the  torque  about  point  W  that  is 
exerted  by  the  mercury  on  its  channel  walls.  Let 

r  -  torque  vector  about  W  of  the  forces  which 
the  channel  exerts  on  the  mercury, 

r.  =  torque  vector  about  W  of  the  centrifugal  (19) 
forces  exerted  on  the  mercury, 

Tp  =  torque  vector  about  W  of  the  Coriolis 
forces  exerted  on  the  mercury; 

then  for  rotational  equilibrium  of  the  mercury  mass  taken  as 
a  whole  we  have 


T  +  tL  +  TC .*  ^IM  ’r*X  VR  da  (20) 

since  W  is  a  fixed  point.  Here  the^jintegral  is  the  moment  of 
momentum,  or  angular  momentum,  and  r  is  a  vector  from  W  to 
the  elementary  mass  dm.  Despite  the  fact  that  the  mercury 
appears  to  be  flowing,  the  moment  of  momentum  remains  con¬ 
stant;  hence  the  right-hand  side  of  Eq.  20  vanishes,  and 
Eq.  20  becomes 


t  =  -  L  ■*  Tp.  (21) 

Noting  Eq.  6  we  see  that  the  centrifugal  force  acting  on 
ar.y  elementary  mass  dm  has  a  line  of  action  which  passes 
through  the  precession  axis,  and  hence  has  a  zero  moment 
about  this  axis.  Integrating  over  the  entire  mass  it  follows 
that 


tL  =  °- 


Finally  it  follows  from  Eq.  8  that 

TC  “  2  i  ^  X  ^XVR^  dm» 


whence 


(22) 


(23) 


Tc  =  -  2  J  J  -  (r"-?)VRl  dm.  (24) 

M 

Since  r  and  £  are  orthogonal  r"0  *  0,  and  the  second  term  in 
the  integrand  of  Eq.  24  vanishes,  leaving 
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TC 


r • V^dm . 


(25) 


Noting  the  figure  below  let  us  place 

r*  =  W<f  +  r*  (26) 

s 

where  WS  is  the  vector  from  W  to  S,  and  r  is  the  vector 

s 


from  S  to  the  elementary  mass  dm.  Substituting  Eq.  26  in 
Eq.  25  we  obtain 

Tc  ■=  -  2/f  f  (WsVs)-VR  dm 

*  -  20  W  *  J  VR  dm  -  2  J  ^  iyVR  dm; 
however  the  first  term  vanishes  due  to  Eq.  16,  leaving 


TC 


dm. 


(27) 


For  either  the  eccentric  annular  shape  or  the  crescent  shape, 
and  for  either  laminar  flow  or  turbulent  flow  the  radial 
component  of  Vr  at  any  element  dm  is  the  negative  of  that 
component  at  the  element  dm  which  is  located  symmetrically 
opposite  that  plane  of  symmetry  of  the  mercury  which  contains 
the  spin  axis,  shown  in  the  following  figure.  It  follows 
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that  the  contribution  of  the  mercury  which  lies  on  one  side 
of  this  plane  to  the  integral  in  Eq.  27  cancels  the  contri¬ 
bution  of  the  mercury  which  lies  on  the  other  side;  hence 
the  integral  vanishes,  and 

Tc  »  0.  (28) 

Equations  22  and  28  in  Eq.  21  now  give 

T  *  0.  (29) 

We  thus  see  that  the  moment  about  the  precession  axis  of  the 
forces  exerted  by  the  channel  upon  the  mercury  is  zero.  It 
follows  that  the  resultant  of  these  forces  passes  through  W, 
and  hence  may  be  considered  to  act  at  W,  as  shown  below. 


to  cepter  of 
gravity  of 
mercury 


F  is  the  force  exerted  by  the  channel  on  the  mercury;  hence 
(-F)  is  the  force  exerted  by  the  mercury  on  the  channel,  and 
thereby  on  the  main  spinning  body,  as  shown  in  the  following 
figure.  (-F)  may  be  replaced  by  the  two  component  vectors 
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gravity 


Fj  and  Fj ,  the  magnitudes  of  which  are,  noting  Eq.  18, 

If  |  *  | F I  sin^  *  M  r  0*  sin*.3*  -  M  y  ft*  , 

1  (30) 

If.  i  *  If  I  cos*,  *  m  r.^a  cos*.*  m  xnP* 

I  10  10 

where  x0  and  y0  are  the  x  and  y  coordinates  of  the  center  of 
gravity  of  the  mercury  measured  with  reference  to  a  Cartesian 
coordinate  system  having  W  as  an  origin  and  WS  as  the  positive 
x  direction.  Let 


WS«{;  (31) 

then  the  moment  of  Ft  about  the  spin  axis  is,  noting  Eq.  30, 


Moment  exerted 
by  the  mercury  m 
about  the  spin  z 
axis 


If,  I  {  -  -  M  y0{/s‘ 

counterclockwise . 


(32) 


This  acta  in  the  same  sense  as  p,  and,  together  with  an  equal 
moment  due  to  the  other  damper,  results  in  a  rate  of  increase 
of  v  given  by 


(33) 


It  should  be  noted  that  y#  is  negative. 

Fj  and  F,  give  rise  to  torque  vector  components  r  and  r  , 
respectively,  about  the  center  of  gravity  of  the  mainxspin-  ” 
ning  body  as  shown  in  the  following  figure,  where  noting  the 
figure  at  the  bottom  of  page  2, 
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W 


S 


y‘ 


=  IFJL, 
-  |F  |L, 


(34) 


L  being  the  distance,  closely,  between  W  and  the  center  of 
gravity  of  the  main  spinning  body.  Since  the  torque  vector 
of  the  external  forces  acting  on  this  body  about  its  center 
of  gravity  is  the  velocity  of  the  tip  of  its  angular  mo¬ 
mentum  vector  H  for  the  center  of  gravity,  we  see  that  the 
torque  vector  components  7X  and  Ty,  when  multiplied  by  2  to 
include  the  effect  of  the  other  damper,  may  be  regarded  as 
velocity  components  of  the  tip  of  this  vector.  Since  the  H 
vector  lies  along  the  precession  axis  and  is  of  magnitude  A py 
closely,  we  see,  noting  Eqs.  30  and  34,  that  T  causes  W  to 
move  toward  S  with  a  speed  A 


Velocity  component  _  _  ;  ^  2  Tx _ 
of  W  toward  S  *  "  Av 

2My0/3*L* 

“  Av  ’ 


(35) 


It  is  this  velocity  that  lines  up  the  precession  and  spin 
axes,  and  hence  gets  rid  of  the  wobble. 

Similarly,  we  see  that  r  causes  W  to  move  in  the  y 
direction  with  a  speed  y 


Velocity  component  _  2  ^  Ty  ^  m  ^ 

of  W  in  y  direction  Av  kv  *  '  * 

the  factor  2  being  put  in  to  include  the  effect  of  the  other 
damper.  This  subtracts  directly  from  the  y  component  of  the 
velocity  of  S  relative  to  W,  the  y  component  of  the  velocity 
of  S  being  due  to  the  precession  and  equal  to  {/J.  The  effect 
of  this  is  to  reduce  the  precessional  angular  velocity  by  the 
amount 
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Reduction  in  precessional  ^  ^  .  . 

angular  velocity  *"  A vl  ^37' 

since  the  angular  velocity  of  S  about  W  is  the  precessional 
angular  velocity. 

Another  effect  of  the  velocity  component  of  W  in  the  y 
direction,  given  by  Eq.  36,  is  to  shift  the  precession  axis 
off  of  the  H  vector  for  the  main  body.  Noting  the  figure 
below  we  see  that  although  the  H  vector  still  passes  through 
W,  the  precession  axis  now  penetrates  the  plane  of  the  damper 


-  V  y* 

ww  = - =*- 

l  v  -  v 
sy 


Here  v  is  given  by  Eq.  36,  and 
wy 


(38) 


WS  -  4,  vsy  -  «/9. 

Attention  has  already  been  called  to  the  fact  that  the 
translational  effects  of  the  two  dampers  cancel. 
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RELATION  OF  THE  VISCOUS  FRICTIONAL  POWER 
LOSS  TO  THE  TORQUE 


The  force  F  exerted  by  the  channel  on  the  mercury  (see 
figure  at  the  bottom  of  page  11)  arises  through  the  action 
of  the  stresses  transmitted  by  the  channel  walls  to  the 
mercury.  These  stresses  consist  of  the  normal  stress,  or 
pressure,  and  the  shear  stress,  which  is  due  to  the  viscosi¬ 
ty  of  the  liquid.  Taking  moments  about  the  spin  axis,  which 
passes  through  S,  we  see  that  since  the  pressure  contributes 
nothing,  the  entire  moment  is  contributed  by  the  shear  stress. 
Noting  Eq.  32  and  the  figure  on  page  12  we  see  that  this  mo¬ 
ment  is  of  magnitude  |F1|{,and,  looking  down  from  the  top, 
acts  clockwise  on  the  mercury  (figure  below).  Here  the  eccen¬ 
tric  annular  shape  is  shown  merely  to  fix  ideas.  At  present 


motion  of  mercury 


we  are  making  no  restriction  upon  the  shape  of  the  mercury. 

Since  the  main  body  rotates  with  an  apparent  angular 
velocity  (p-v)  clockwise,  it  delivers  a  power 

P  -  IF,  l{  (A-O  -  -  M  y9(fl*  Ot-P)  (39) 


to  the  mercury.  Although  the  mercury  is  apparently  flowing, 
its  kinetic  energy  remains  constant;  the  potential  energy 
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which  it  has  by  virtue  of  the  potential  function,  Eq.  7, 
which  arose  because  of  the  centrifugal  force,  is  constant; 
and  the  power  delivered  by  the  Coriolis  force  field  is  zero, 
since  the  Coriolis  force  is  everywhere  perpendicular  to  the 
velocity  of  the  mercury.*  It  follows  that  the  entire  power 
(Eg.  39)  is  dissipated  as  heat,  and  is  equal  to  the  viscous 
frictional  power  loss  in  the  mercury I  Substituting  feq.  39 
in  Eq.  35  we  obtain 


Velocity  component  ;  2PL* 
of  W  toward  S  »  -  5  *  Hp£(fi-v)  * 


(40) 


hence  the  viscous  friction  loss  P  is  the  primary  factor  in 

the  determination  of  k  »  and  hence  in  the  elimination  of  the 
wobble. 


2  Instead  of  using  the  constancy  of  the  potential  energy 
to  show  that  the  power  delivered  by  the  centrifugal  force 
field  to  the  mercury  is  zero,  we  could  also  compute  this 
power  directly;  thus  noting  Eqs.  6  and  7  we  have 

Power  delivered  to  r  r 

mercury  by  centrifugal  ■  -  \  pVR*V*dv  «  p  \  [*V*VR-V*  (#VR) ]dv 

force  field  Jv  Jv 

where  p  is  the  density,  and  the  integral  extends  throughout 
the  volume  of  the  mercury.  Applying  the  divergence  theorem, 
and  noting  that  V*VR  ■  0  since  the  fluid  is  Incompressible, 
this  becomes 

Power  delivered  to  p  , 

mercury  by  centrifugal  •  -  p  \  #V«  *  da  *  0 

force  field  Ja 

where  3aT  is  the  vector  element  of  area,  and  the  integral  ex¬ 
tends  over  the  mercury  surface;  for  Vr*3T  -  0  because  at  this 
surface  VR  is  tangent  to  the  surface,  and  hence  perpendicular 
to  the  normal. 

The  fact  that  in  the  case  of  the  eccentric  annular  shape 
the  mercury  surface  is  not  simply  connected  causes  no  diffi¬ 
culty  in  applying  the  divergence  theorem,  for  a  plane  cross¬ 
cut  could  be  introduced  to  render  the  surface  simply  connected 
without  altering  the  values  of  the  integrals. 
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MINIMUM  PRINCIPLE  FOR  THE  VISCOUS  FRICTIONAL 

POWER  LOSS 


Since  the  mercury  is  incompressible  =  0,  and  at  any 

instant  the  lines  of  flow  cannot  terminate  anywheza,  but  must 
be  closed.  It  follows  that  at  any  instant  the  flow  may  be 
considered  to  be  composed  of  an  infinite  number  of  closed 
elementary  tubes  of  flow,  figure  below.  Let  us  suppose  that 


the  flow  is  stationary — that  is,  the  flow  pattern  does  not 
change  with  time,  it  being  remembered  that  while  observing 
this  pattern  we  are  rotating  counterclockwise  with  an  angular 
velocity  0.  We  shall  now  show  that  if  the  flow  is  stationary. 
and  if  the  shape  of  the  velocity  field  is  held  fixed,  so  that 
the  lines  and  tubes  of  flow  are  unaltered,  then  the  flow  will 
distribute  itself  among  the  various  tubes  in  such  a  wav  that 
the  viscous  frictional  power  loss  is  a  mlnimumT 

Let  u_  choose  a  set  of  rectangular  xyz^coordinate  axes 
to  measure  the  apparent  motion  of  the  mercury;  and  let  vx, 

Vy,  and  v,  denote  the  x,  y,  and  t  components  of  Vr,  respec¬ 
tively;  then  the  normal  and  shear  stresses  in  the  liquid  which 
arise  due  to  viscosity  are,  respectively,* 


Rouse,  Hunter.  Fluid  Mechanics  for  Hydraulic  Engi¬ 
neers.  New  York,  McGraw-Hill,  1938.  Pp.  144*45. 


1? 


(41) 


The  corresponding  time  rates  of  change  of  strain  are,  re¬ 
spectively, 


€  =  - C  £  =  - JL  p  * 

x  8x  }  y  8y  *  z 


B^z 

8z  * 


8v 

8  v.. 

»  • 

y  —  y  - 

xy  'yx 

x 

8y 

y 

8x  ’ 

8v„ 

8v 

•  • 

V  «  V  « 

rzx  xz 

z 

8x 

+ 

X 

8c  * 

(42) 


Noting  that  the  fluid  is  incompressible  the  viscous  fric¬ 
tional  powerless  per  unit  volume  is  hence 


r  /0v.  8v„\* 


L\  8y  Ox 


Power  loss  per  +(??£)*  +(!!£)* 

unit  volume  Ox/  \  By  J  \  9z  J 

(£•*£)'] 


(43) 


Integrating  throughout  the  volume  of  mercury  we  obtain  for  P, 
the  total  viscous  frictional  power  loss 


M  f  iJ/®vxV  .  /®vyV  .  /•v«\,l 

“  jv  UW  +  {~WJ  +  \yr)  J 


(44) 


Holding  the  shape  of  the  velocity  field  fixed  let  us  give 
VR  a  variation  5VR,  the  corresponding  variations  in  vx,  v  , 
and  v,  being  $vK,  6vy,  and  dvz,  respectively.  At  every  point 
«VRtfi  us  has  the  same  direction  as  hence  it  is  really 
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|VR|  which  is  being  changed.  We  require  also  that  the  varied 
field  satisfy  the  actual  surface  boundary  condition  that  at 
the  solid  boundary  surface  provided  by  the  channel  the  fluid 
velocity  be  that  of  the  solid  surface.  We  thus  have  at  the 
mercury  surface 


Vr  =  relative  velocity  of  solid  bounding 
surface  at  such  a  surface , 

6Vr  =  0  at  a  solid  bounding  surface,  (45) 

Shear  stress  tangent  to  a  free  surface  =  0 
for  the  actual  velocity  field. 


The  variation  in  P  can  now  be  obtained  from  Eq.  44,  thus 


6  P 


=  *  1  {< 


r8vx  *  8vx  + 

aT  4  W  + 


8y  8  y  8  z 


8vzl 

T5FJ 


/®V  9vy 
+  (iy  + 


8v 

6  -K-—  +  6 

By 


8vy\ 

S 


/0v, 


+ 

V  8z 


9v 


7. 


+  - 


\ 

/®V. 


0y 

8v 


8v  0v  s 

4  "8#  +  5  77 


V  8x 


590 


8v 


8v 


4  -sr  +  4  -sf. 


dv 


lince 


8vx  »  8vv 

*  "5x‘  ”  57  8vx>  8  W 


y  *y 

it  follows  that  Eq.  46  may  be  written 

8v_  8v 

u\ 

v 

/< 


8  , 
87  6v 


etc. 


«p-  2,1  ^  ^  «vy  +  -5F  H  4vz] 

(^7  +  -St)  (4;  4vx  +  fx  Svy) 

+  T?)  (n  4vy  +  ^  4vz) 


V  8  x 


+  (-*  +  (£  *VZ  +  £  4vx)|  dv. 


(46) 


(47) 


(48) 
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Considering  only  the  first  term  in  the  integrand  of  this 
integral  we  have 


f  f  f  0VX  a 

4"  j  j  j  87  6vx  dydzdx- 


(49) 


y  z  x 

Integrating  by  parts  this  becomes 


■{  ttf 


*,x*  (y,z) 


•— £  6v  I 
8x  xjx1 vy,z> 


8  ^ 

dydz  -  ^  5vx  dydzdx^  (50) 


z4  8x*  x 


where  xt(y,z)  and  x*(y,z)  are  values  of  x  on  the  surface. 
There  may  be  more  than  two  such  values,  but  this  causes  no 
difficulty.  Denoting  the  direction  cosines  of  the  normal  to 
the  surface  by  coso^ ,  cos.8  ,  cosy  Eq.  50  becomes 


sf  5vx 


f*  &  v  'I 

cosot  da  -  \  - ~  5v  dvr  (51; 

1  Jv  8x  x  J 


where  the  first  integral  extends  over  the  surface  of  the 
meicury,  and  the  second  extends  throughout  its  volume.  Pro¬ 
ceeding  in  the  same  manner  with  the  other  terms  in  the  inte¬ 
grand  of  Eq.  48,  Eq.  48  becomes 


<5?  “  2#Jf 


f  f  r8v*  8vy  ®va 

*  j  \ 2  Sx  SVX  cosai +  w  6vy  cos^>  4v*  cosr» 

d  *■ 

8Vy\ 

+  VSy"  +  Try  (Svx  C0S<,1  +  Svy  C0SV 


>vv  »v„ 


+  [nr  +  WJ  (6vy  COSri  +  dvz  cosV 

+  (~w  *  tt)  (Svz  ccsai  +  Svx  cosri)}  da 

f  f  r*,y*  0,vy  ®*vz 

-  2u  V  -  2 - *  «v  +> — J  *v  + - -  «v 

^v  ^  8x‘  x  ay*  y  a** 

,  /  8  x  8  'N  /*vx  .  8vy\ 

+  (dvx  8?  +  4vy  n;  V?r  W 


(52) 
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8 

8z 


+  5v 

z 


+  !v\ 

l  82  8  yj 


+  (  6v 


8 

8x 


8v. 


+  Sv 


jl\  f: 

8zJ  l,  8x 


+ 

+  TT 


)} 


dv. 


In  view  of  Eq.  41,  Eq.  52  may  be  written 


<5P 


2  j*  ^°x  cos<*i  +  TXy  cos^l  +  Txz  cosyl)«v) 


+  (r  cosa  +  cr  cosp,  +  r  „  cosy.)6v 
xy  l  y  l  yz  l  y 

+  (t  cosa  +  r  cos#,  +  <r  cosy,  )5v  Ida 
xz  i  yz  l  z  '  l  z 

V  v  u  fix1 


(53) 


[>  Gf 

<; 


8y 


8y 

0*V 


+  -*VjL 

®y /  ®x  V®y 

9 


y-+lf5K 


+  L2  ("57  +  fr)+*?  (t^+T5t).4v*}  dv- 


But  the  coefficients  of  $vx>  6vy,  and  5vz  in  the  surface 
integral  are  the  components  of  surface  stress;  and  in  the 
volume  integral  we  note  that  since  the  fluid  is  incompres¬ 
sible 


whence 


(54) 


(55) 


21 


NAVWEPS  REPORT  8611 


It  follows 


6  P  =  2 


+ 


+ 


where  S  is 
wr itten 


that  Eq.  53  may  be  written 


the  surface  stress  vector. 


(56) 


Equation  56  may  be 


<5  P 


S*<5VRda 


2 P 


V'6V  Rdv* 


(57) 


But  5Vr  vanishes  at  the  channel  walls,  and  the  tangential 
component  of  S  vanishes  on  the  free  surface,  which  causes 
S*6Vr  to  vanish  there  also.  It  follows  that  the  integrand 
of  the  surface  integral  vanishes  over  the  entire  mercury 
surface;  hence  this  integral  vanishes,  and  Eq.  57  reduces  to 


6P 


2 P 


VR),<5VRdv- 


(58) 


The  condition  for  dynamic  equilibrium  in  the  liquid  is 
PVWR  =  -  Vp  -  pV*  -  2 p/f  x  VR  +  pV*VR  (59) 


where  p  is  the  density  of  mercury,  and  p  is  the  pressure. 

Here  the  left-hand  side  is  the  mass  times  acceleration  per 
unit  volume,  and  the  terms  on  the  right-hand  side  are  the 
force  contributions  per  unit  volume  due  to  the  fluid  pres¬ 
sure,  the  centrifugal  force  field  (Eq.  7),  the  Coriolis  force 
field  (Eq.  8),  and  the  viscous  frictional  force,  respectively. 
Multiplying  Eq.  59  by  and  integrating  throughout  the 
volume  of  the  mercury  we  obtain 


P 


vR)-<5VRdv 


(pVR*WR  +  Vp  +  pV* 


+2  p’fi  *VR)-dVRdv. 


(60) 
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In  order  to  evaluate  this  integral  let  us  consider  the  volume 
to  be  broken  up  into  elementary  tubes  of  flow,  corresponding 
to  which  we  let  the  element  of  volume  dv  be 

dv  =  dtrds  (61) 

where  ds  is  the  length  and  dtr  is  the  cross-section  of  an 
elementary  tube,  as  shown  in  the  following  figure.  Also  let 
the  volume  of  mercury  which  flows  across  any  cross-section 


of  flow 


of  the  tube  per  unit  time  in  the  direction  of  the  unit  tan¬ 
gential  vector  ls ,  in  tfe  flow  5VR  be  d6q;  then 


(62) 


Multiplying  by  Eq.  61  now  gives 

SVrjdv  =  1  dsdfiq, 

K  S 

whence  Eq.  60  becomes 

P  j  (V^.SVj^dv  »  dfiq  ^  ^1S*(VR,WR) 


(63) 


all  line  of 
tubes  flow 


+  ls*Vp  +  pl$*V*  +  2pi s'J  x  VR]ds . 


(64) 


But  the  triple  product  1-’/*  x  VR  vanishes  since  lg  and  VR 
have  the  same  direction; salso 


s 


1  -Vpds 


0, 


(65) 
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s 

Substituting 


p ls*V$ds  “  j  P  ^7  ds  =  P 
in  Eq.  64  and  noting  Eq. 


=  0. 

we  obtain 


SP  =  -  2  J  d<5q  £  Pls*(VR*WR)ds. 

all  line  of 
tubes  flow 


But 


VWR 

VWR 


'VVWR 


(66) 


(67) 


(68^ 


where  R  is  the  radius  of  curvature  of  the  line  of  flow,  and 
1«  is  a  unit  principal  normal  vector  pointing  toward  the 
center  of  curvature  of  this  line.  Substituting  in  Eq.  67 
now  gives 


«  -  -  2  )  J  sr- ds 

all  line  of 
tubes  flow 

=  -  P  y  ddq  £  d (  |  VR I  * )  ,  <5P  =  0.  (69) 

all  line  of 
tubes  flow 


The  actual  value  of  P  is  therefore  stationary  when  the  veloci¬ 
ty  field  is  given  the  variation  dVp. 

In  order  to  show  that  this  value  of  P  is  actually  a 
minimum  we  consider  the  second-order  term  which  was  not 
included  in  Eq.  46,  namely 
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This  term  is  evidently  positive,  and  is  the  viscous  fric¬ 
tional  power  loss  which  would  be  obtained  with  the  velocity 
distribution  6V R  alone.  Since  the  first-order  terms  vanish, 
as  indicated  by  Eq.  69,  the  increase  in  P  due  to  6V R  is 
given  by  the  positive  quantity  (Eq.  70).  Since  this  vanishes 
if  and  only  if  oVr  =  0,  we  see  that  P  takes  its  smallest 
value  when  dVR  =  o,  and  hence  can  only  be  increased  by  devi¬ 
ating  from  the  actual  velocity  distribution  V„  in  the  manner 
described. 

As  an  application  of  this  principle  of  least  power  loss 
we  see  immediately  that  in  the  case  of  the  eccentric  annular 
shape  the  relative,  or  apparent,  velocity  of  the  mercury 
will  in  some  places  be  greater  than,  and  in  other  places 
smaller  than  that  of  the  channel.  It  follows  that  an  ob¬ 
server  riding  with  the  main  body  would  see  the  mercury  going 
in  the  direction  of  motion  of  the  channel  in  some  places, 
and  in  the  opposite  direction  in  others.  This  is  the  situa¬ 
tion  which  pertains  in  a  wave  motion. 

As  another  application  we  see  that  since  P  is  a  minimum 
for  the  actual  distribution  of  VR  it  is  not  sharply  dependent 
upon  this  distribution.  It  follows  that  if  an  approximate 
distribution  of  VR  is  obtained  by  some  approximation  proce¬ 
dure,  the  accuracy  obtained  for  P  is  much  greater  than  that 
obtained  for  VR.  This  is  important,  for  the  quantity  in 
which  we  are  really  interested  is  P,  since,  as  we  saw  in 
Eq.  40,  it  is  P  which  determines  the  rate  of  damping  of  the 
wobble .  Finally  we  note  that  the  approximate  value  of  P 
obtained  is  slightly  too  large,  since  it  would  be  reduced  if 
the  actual  distribution  of  VR  instead  of  the  approximate  one 
were  used  in  computing  it. 


IMPROVED  TREATMENT  OF  THE  ACTION  OF  THE  DAMPER 
ON  THE  MAIN  SPINNING  BODY 


The  effects  of  the  forces  exerted  by  the  mercury  on  the 
channel,  and  hence  on  the  main  spinning  body,  were  previously 
considered  (pp. 13-14)  with  regard  to  their  effect  in  modify¬ 
ing  the  motion  of  this  body.  The  treatment  was,  of  course, 
approximate;  nevertheless  it  was  adequate  to  show  how  the 
damper  reduces  the  precession,  or  wobble,  how  it  reduces  the 
precessional  angular  velocity,  and  how  it  alters  the  preces¬ 
sion  axis.  It  is  open  to  tne  objection,  however,  that  the 
forces  exerted  by  the  mercury  were  computed  using  the  motion 
which  the  body  had  prior  to  being  acted  upon  by  these  forces. 
Although  this  situation  suggests  a  method  of  successive 
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corrections  wherein  we  first  recompute  the  forces  using  W 
instead  of  W (see  figure  on  page  14),  and  the  new  value  of1 
precessional  angular  velocity;  and  then  compute  second  cor¬ 
rections  to  the  motion  of  the  main  body--and  so  on;  nevez- 
theless  we  shall  avoid  this  by  taking  immediate  account  of 
the  fact  that  with  the  damper  acting  the  angular  momentum 
vector  H  of  the  main  spinning  body  does  not  lie  along  the 
precession  axis. 

Noting  the  figure  below  and  the  figure  on  the  bottom  of 
page  2  we  see  that  for  W5  =  |  small  the  K  vector  is,  closely, 


composed  of  a  component  of  magnitude  hv  lying  along  the  spin 
axis,  and  a  component  of  magnitude  at  right  angles  to 
this  axis  (figure  below).  Denoting  the  point  where  the  line 


m 

L 


extending  from  the  center  of  gravity  of  the  main  body  along 
the  H  vector  penetrates  the  plane  of  the  damper  by  H,  which 
should  cause  no  ambiguity,  we  see  that 
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HS  = 


UL S 

__k_  l  =  Ml 

A  y  D  A>/  * 


(71) 


Since  the  velocity  of  the  tip  of  the  H  vector  is  the 
due  to  the  external  forces,  in  this  case  those  exerted  by 
the  mercury,  it  follows,  noting  the  figure  on  page  13,  the 
fact  that 


WH  -  £  (l  -  (72) 

and  the  fact  that  tx  and  t y  must  be  multiplied  by  2  in  order 
to  include  the  effect  of  the  other  damper,  that 

2  17,1  *)--*£««.  03) 

2  |Tyl  "  «  C1  -  S?)  f  (t)-¥  (A,'-^)-  (74) 

In  connection  with  Eq.  73  we  note  that  although  W  is  not 
stationary,  it  is,  nevertheless,  the  instantaneous  center  of 
S  in  its  motion.  The  velocity  component  of  S  in  the  direc¬ 
tion  W5  is  therefore  zero;  hence  the  component  of  the  velocity 
of  point  H  in  this  direction  is  the  same  as  the  component  of 
the  velocity  of  point  H  relative  to  S  in  this  direction. 
Equation  73  therefore  follows  from  Eq.  71. 

If  we  let  A  denote  the  decrease  in  precessiona.l  angular 
velocity  due  to  the  dampers,  thus 

A  =  —■  -  /9,  (75) 

then  Eq.  74  becomes 

2  I7y|  -  ^£5.  (76) 

Substituting  Eqs.  30  and  34  in  Eqs.  73,  74,  and  76,  we 
obtain,  respectively, 

~  ({/>>  -  2  (*£)  A„,  07) 

{(Av-B0)  -  2ML*0xo,  (78) 
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£AB  2ML,0xq. 


(79) 


Further  consideration  of  these  equations  must  await  the 

determination  of  x„  and  v  . 

o  }  o 


ECCENTRIC  ANNULAR  CONFIGURATION- -NATURE  OF  THE 
FREE  SURFACE  OF  THE  MERCURY.  SPECIFICATION 
OF  THE  SHAPE  0?  THE  VELOCITY  FIELD 


We  saw  earlier  that  a  given  mercury  damper  may  operate 
in  several  configurations  (pp.  2-7)  depending  upon  the  size 
of  the  wobble,  quantity  of  mercury,  md  dimensions  of  the 
channel;  also  we  saw  that  one  mercury  damper  may  not  exper¬ 
ience  the  same  sequence  of  configurations  as  <•  not  her.  How¬ 
ever,  in  all  cases  the  final  configura  ico  i0  the  same  for 
all  mercury  dampers,  namely,  eccentric  annular.  The  deter¬ 
mination  of  the  nature  of  the  final  damping  of  the  wobble 
hence  in  all  cases  requires  the  analysis  of  the  eccentric 
annular  configuration.  We  must,  in  particular,  determine 
the  viscous  frictional  loss  P,  nee  we  saw  from  Eq.  40  that 
it  is  this  quantity  which  is  of  primary  importance  in  damping 
out  the  wobble.  In  determining  P  we  shall  apply  the  minimum 
principle  derived  previously  (pp  17-25)  to  determine  the 
velocity  distribution;  however,  in  order  to  oo  this  we  must 
f^rst  specify  the  shape  of  th*  velocity  field. 

We  saw  on  page  5  that  if  there  were  no  wall  friction  the 
inn._r  free  mercury  surface  would  be  a  circular  cylinder 
centered  at  W  as  shown  in  the  figure  on  page  5.  With  wall 
friction,  however,  this  inner  surface  is  no  longer  exactly 
a  circular  cylinder,  and  is  no  longer  centered  on  W,  but  is 
displaced  as  indicated  in  the  figure  on  page  15.  Since  for 
cur  purpose  the  circular  shape  offers  no  particular  advan¬ 
tage,  this  deviation  from  it  is  no  cause  for  concern. 

Using  polar  coordinates  r  and  9  let  us  consider  the 
function 


r  «  a  +  b  cos 9  +  c  sine  (80) 

where  b  and  c  are  small  compared  with  a,  as  shown  in  the 
figure  on  the  following  page.  Equation  80  may' be  written 

r*  *  ar  +  bx  +  cy;  (81) 
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whence 


r*  -  2ar  +  a1  =  -  ar  4-  a*  4  bx  +  cy, 
(r-a)a  =  -  ar  t  a*  +  bx  +  cv. 


(82) 


Retaining  only  first-order  terms  in  b  and  c  we  see  from 
Eq.  PO  that  (r-a)*  is  of  second  order  in  these  quantities; 
hence  the  left-hand  *ide  of  Eq.  82  may  be  replaced  by  zero, 
giving 


r 


(83) 


Squaring  and  omitting  terms  of  order  higher  than  first  ordei 
in  b  and  c,  we  obtain 

x*  +  y*  *  a*  +  2bx  4-  2cy; 

or,  again  good  through  first-order  terms, 

(x-b)1  +  (y-c)*  -  aa .  (84) 

We  thus  see  that  good  through  first-order  terms  in  b  and  c 
the  polar  cu^ve  (Eg.  80)  is^a  circle  of  radius  a  whose  center 
is  at  point  x  *  b .  v  »  c .  to  this  accuracy  Kq.  80  may  "Ee" 
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used  in  place  of  the  equation  for  a  circle,  and  for  our 
purpose  is  much  more  convenient. 

For  the  present  we  shall  confine  our  attention  to  the 
jase  where  the  eccentricity  (distance  of  center  from  S)  of 
the  inner  mercury  surface  is  small,  which  corresponds  to  tne 
final  phase  in  the  operation  of  any  mercury  damper.  It  is 
evident  that  any  particle  of  mercury  which  is  in  contact 
with  the  channel  wall  moves  in  a  circle  of  zero  eccentricity, 
since  it  must  have  the  same  velocity  as  that  point  on  the 
wall  with  which  it  is  in  contact.  Supposing  as  an  approxi¬ 
mation  that  all  of  the  particles  of  mercury  move  in  circular 
paths  it  follows  that  the  eccentricities  of  these  paths  in¬ 
crease  from  zero  as  we  move  away  from  the  channel  walls  and 
into  the  mercury.  In  particular  we  see  that  the  central 
part  of  the  free  surface  of  the  mercury  bulges  outward  so  as 
to  Increase  the  depth  of  the  mercury  where  it  is  deepest 
(radially),  and  shrinks  inward  so  as  to  decrease  this  depth 
where  it  is  most  shallow.  The  free  surface  is  hence  by  no 
means  cylindrical.  Keeping  these  facts  in  mind  and  noting 
Eq.  80  we  shall,  as  an  approximation,  choose  the  function 


r 


X (ij -a)  —  X  b 


(85) 


to  specify  the  shape  of  the  velocity  field.  Here  cylindrical 
coordinates  are  used,  the  origin  being  at  S,  the  center  of 
the  channel  (as  shown  in  figure  on  the  following  page);  also 

r  =-  outer  radius  of  channel. 

2 

z  "  half  the  axial  (vertical)  width  of  the 
A  channel , 

(86) 

a  =  radius  of  inner  free  surface  with  no  v  ' 

wobble , 

b  =  worst  eccentricity, 

X  =  parameter. 

For  any  values  of  X,  b,  and  z  Eq.  85  gives,  accurate  through 
first-order  terms  in  b,  a  horizontal  circle.  As  the  para¬ 
meter  X  varies  from  zero  to  one  with  b  and  z  fixed  a  family 
of  circles  is  obtained,  the  first  being  on  the  outer  channel 
wall,  and  the  last  lying  cn  the  free  surface.  When  z  *  *  zJf 
or  X  -  0,  corresponding  to  channel  walls,  the  eccentricity  of 
these  circles  is  zero.  When  X  ■  1  Eq.  85  becomes  the  equation 
of  the  free  surface,  shown  in  the  figure  on  the  following  page. 

The  circle  which  lies  at  the  center  of  the  free  surface,  for 
which  z  «  0,  X  ■  1,  has  the  greatest  eccentricity,  namely  b. 
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r  -a 


The  total  volume  of  mercury  is 


Volume  of 
mercury 


(87) 


31 


\ 

% 


NAVWEFS  REPORT  8611 


Volume  ot„ 
mercury 


-  *  a  4  va 

2rzi  [Va  -  15  b  j- 


(88) 


Accurate  through  first-order  terms  in  b  we  thus  hive 

Volume  of  a  2 
mercury  l v  a  ' 


(89) 


From  Eq.  89  we  see  that  accurate  through  first' order  terms 
in  b  the  radius  a  of  the  circles  which  compose  the  free  sur¬ 
face  of  the  mercury  is  the  same  as  that  of  the  inner  free 
cylindrical  surface  of  the  mercury  when  there  is  no  wobble. 

Referring  to  the  figure  on  the  previous  page  we  see  that 
the  y  and  z  coordinates  of  the  center  of  gravity  of  the  mer¬ 
cury  are  both  zero  due  to  symmetry.  The  x  coordinate  x  times 
the  mercury  volume  is 


x  times 
mercury 
volume 


-  3ab 


V*1J  , 

- }  1  cos  S 
J  J 


I  1 1  It.  f  cos9  dzd9dr 

•Zj  o  a-b]  x-  'cosfl 

ifM  "  ffYlcos* 

3  J-zx  Jo  1  *  L  \*l/  j 

-  3ab3  [l  -  ( —)  1  cos*  9 

L  \3:y  J 

+  b8  ["  1  -  (•£-)  1  cos8??*  co$8dzd8 

L  \zkJ  J  J 

f  j  ;■  ft  t>  -  (■&)  *  i  -■[>  • = (0 

*  ■  (0  -  0}  )*• 

o.K,  r,*  A  4  b*  ?  4.  3  lY 

2i  bz  a  1  »  r  t  ~r*  i  1  -  r  t  r  -  r 

11  V  V  4  \  35  VJ 


(90) 


x  times  A  *  4  a\ 

mercury  *  IrbCj  f ~  a  +  b  j. 


(91) 


volume 


Dividing  by  the  mercury  volume  (Eq.  88)  we  now  obtain 
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2b(35a14-6b1) 

7  [15 (r* -a*  )-4b*  ] 


(92) 


Accurate  through  first-order  terms  in  b  we  thus  have 


x  « 


2ha 


3>ra-a  ) 


(93) 


The  radius  of  curvature  of  the  curve 

r  «  f(tf) 

in  polar  coordinates  is 


(94) 


R 


(r*+r ?1  )a 
rS+2r  f,-rr" 


(95) 


where  a  prime  indicates  differentiation  with  respect  to  ©. 
Applying  this  to  the  function 


r  *  a„  +  cos © 
o  o 


(96) 


where  aQ  and  b0  are  constants  we  obtain 


(a*  4  2a.  bA  cos ©  4-  b1  cos1©  4-  b*  sin*©)1 
_ 0  0  0  0  0 

a* 4-2*  b  oos©4-b*cos* ©4-2b*  s in* ©4-a  b  cos©+b*  cos* © 
oco  o  o  oo  o 


-1  b^  !  co»»)! 

4  3aobo  c05*  4  2bJ 


(97) 


Applying  the  binomial  theorem  now  gives 
s  x 

R  -  4|  (bj+2*#bec05») 


4il  <•!>**  2  <»:«•.», co»»>*  +  ...J 
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x  [(aj)”1  +  (-1)  (aj)"a  (3a0bQcose+2bJ) 

+  (-l)(-2)(aj)“3  \  (3aobocos0-f2ba)2  +  ...] 


R  =  a0  +  b 


+  b 


+  t  ao 


R  =  a0+^r 
0  a0  u 


0  (a  3aJ  cos0  -  ao  -T  3a0  cos0 

ao  ao  J 

o  (-4+  <-»<-«  -«  i9a»  cos‘e 

0  0 

(-1) 

'  a  ' 

o 

+  !ifi4*!  cos’e)  4 


+  2  a02a0  cos» 


3a,,  cos5 
o 


2 

✓  a  J 

o 


3 


-  2  +  ^  (  9-9+  -|J  cos 2 6 


i) 


+  . . . 


(98) 


R  =  a.  +  *s4—  (3cos*0-l)  +  ...  . 
o  2a0 

Accurate  through  first-order  terms  in  b  we  thus  have 


R  «  a0. 


(99) 


Comparing  Eqs,  85  and  96  we  can  apply  this  result  to  Eq. 
85,  and  thereby  obtain 


R  “  rf  -  A(rJ-a) 

accurate  through  first-order  terms  in  b. 


(100) 


ECCENTRIC  ANNULAR  CONFIGURATION- -MAGNITUDE 
OF  THE  VELOCITY  FIELD 


Corresponding  to  dAds  we  obtain  a  tube  of  flow  whose 
cross  section  is  an  elementary  parallelogram.  In  this  tube 
the  ranges  of  values  of  A  and  z  are  infinitesimal;  however, 

0  vanes  from  zero  to  2»,  and  the  shape  of  the  parallelogram 
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tube 


varies  with  0,  as  shown  in  the  above  figure.  S’nce  the  shape 
of  the  velocity  field  is  completely  specified  by  Eq.  85;  the 
volume  of  mercury  which  flows  in  any  tube  dAdz  per  second  can 
be  specified  by  the  velocity  function  q(A,z),  th”.s 


Volume  per  second 
in  tube  dAdz 


q(A ,z)dAdz . 


(101) 


Since  this  volume  per  second  is  the  same  in  all  parts  of  the 
tubet  the  velocity  at  any  point  can  be  obtained  by  dividing 
Eq.  101  by  the  tube  cross  section.  Accordingly  we  shall  now 
determine  this  cross  section. 


The  altitude  of  the  elementary  parallelogram  which  com¬ 
poses  this  cross  section  is  dz.  Writing  Eq.  85  in  the  form 

r  -  f(0,A, s)  (102) 

we  see,  looking  down  from  the  top,  that  the  base  of  the 
parallelogram  is  the  distance  between  the  curves  (Eq.  102) 
corresponding  to  A,z  and  A  +  dA,  s,  respectively.  This 
distance  depends  upon  0,  and  is  measured  along  an  orthogonal 
trajectory  of  the  curves,  as  indicated  in  the  schematic 
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diagram,  above.  The  values  of  dr  and  d8  corresponding  to 
this  distance  are  given  by  the  relations 


dr 


8r 

To 


da  + 


(103) 


1  dr  r  _ 

r  d$  br  * 

89 


(104) 


in  which  r  is  given  by  Eq.  102.  Here  Eq.  103  is  obtained  by 
differentiating  Eq.  102,  and  Eq.  104  is  the  condition  for 
orthogonality.  These  equations  may  be  written 


.  8r  ..  8r  .. 
dr  -  57  da  -  dx, 

dr  +  r* da  ■  0, 


(105) 


which  constitute  a  set  of  two  linear,  algebraic  equations  in 
dr  and  da.  Solving,  we  obtain 
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2  8r 


(106) 


It  follows  that  ds  ,  the  base  of  the  parallelogram,  is  given 
by  11 


ds 


n 


/  (dr)’+(rd9)’ 


r 

0r 

8A 

'I8M 

y7r7W 


dA. 


ds 


n 


since  we  see  from  Eq. 


85  that 


dr 

8X 


is  negative. 


(107) 


Dividing  Eq.  101  by  dsndz,  the  cross  section  of  the  tube, 
using  Eq.  107,  we  obtain  for  the  velocity  of  the  mercury  in 
the  tube 


V 


s 


dr 

VX 


(108) 


Here  r  is  given  by  Eq.  85,  and  q(\,z)  will  later  be  deter¬ 
mined  so  as  to  minimize  the  viscous  frictional  power  loss  P. 

Expanding  the  radical  in  Eq.  108  by  the  binomial  theorem 
gives 

ex 

Substituting  fro*  Eq.  85  and  retaining  only  terms  through 
first  order  in  b  we  obtain  finally 


!**•?*.  ***>■: 
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Vs  * 


1  -  ( -Malcos0 

r  -a+2Ab 

a 

L  vv  j 

A*b[l  -  f-^-)  sin® 

1  +  if. _ L _ _ 

2  \  r  ✓  \i 

ra“^^rs“a^“^  ^  1"(T~j 
3  3  L  \ZV 


V  a  aiioAl  J'i  .  _»b  [x  .  (^Vlcossj. 

s  ra~a  l  ra-a  |_  \*x)  J  J 


COS0 


(110) 


CALCULATION  OF  THE  VISCOUS  FRICTIONAL  POWER  LOSS  P 
FOR  GIVEN  b  AND  q(A,z) 


Referring  tc  Eq.  43  we  see  that  in  order  to  obtain  the 
viscous  frictional  power  loss  per  unit  volume  we  must  deter¬ 
mine  the  partial  derivatives  of  each  of  the  velocity  compo¬ 
nents  Vx,  Vy7  and  Vz  with  respect  to  x,  y,  and  z,  respectively. 
We  may,  however,  place  the  coordinate  axes  in  any  way  we  choose. 
Let  us  place  these  with  the  origin  at  the  point  at  which  we 
wish  to  determine  the  loss  per  unit  volume,  the  z  axis  vertical, 
and  the  x  axis  tangent  to  the  line  of  flow,  the  positive 
direction  being  that  of  the  flow.  The  positive  y  axis  then 
extends  outward  normal  to  the  line  of  flow,  as  shown  below. 
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With  this  orientation  of  axes  Vz  =  0 ;  hence 


av 


8V_ 


=  0, 


* 


av. 


=  o. 


8x  ’  8  y  8z 

Furthermore,  since  the  fluid  is  incompressible 


(111) 


av  avv  av 

v'v  -  -5T  +  V  +  T#  =  0i  (112) 

hence  since  the  last  term  vanishes  due  to  Eq.  Ill  it  follows 
that 


8VV  ^  Vy 

8y  “  "  3x'* 

It  thus  remains  for  us  to  compute 

av  av  av  avu  avu 

v  x  X  y  y 

Tx 1  TF>  ~a x  ’  TF* 

We  shall  now  consider  these  separately. 


(113) 


av 

Determination  of  Let  s  be  the  distance  measured 

from  the  origin  along  a  line  of  flow,  and  let  a  denote  the 
angle  which  any  arbi.tr ary  line  of  flow  makes  with  the  posi¬ 
tive  x  direction,  as  shown  in  the  figure  below.  Noting  that 

V  is  a  function  of  position  we  see  that  since 
s 


V  *  V  cosa 
x  c 


(114) 


39 
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it  follow* s  that 


8V 
_ x 

8x 


8V 

~Bx  c0$a 


sina 


da 

8x* 


(115) 


At  the  origin,  which  is  the  only  point  :in  which  we  are 
interested,  a  «  0;  and  Eq.  115  becomes 


8V  8V 

«  _ ? 

8x  8x  * 


(116) 


Since  this  is  merely  the  directional  derivative  of  Vs  in  the 
x  direction,  which  at  the  origin  is  the  direction  of  increas¬ 
ing  s,  it  follows  that 


8V 

s 

~5T' 


(117) 


9Vy 

Determination  of  Referring  to  the  figure  on  page  39 

we  now  have 


V  *=  -  V  sina 

y  s 


(118) 


8x 


sina 


cosa 


ha 

8x* 


At  the  origin  a  =  0,  and  this  becomes 


(119) 


da 

8x‘ 


(120) 


But  ^  is  the  directional  derivative  of  a  in  the  x  direction, 

which  at  the  origin  coincides  with  the  direction  of  increas¬ 
ing  s;  hence  at  the  origin 


8a  8a 

H  71 


curvature  of  line  of  flow. 


(121) 


Equation  120  thus  becomes  finally 


8x 


(122) 


where  R  is  the  radius  of  curvature  of  the  line  of  flow. 
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8V, 

Determination  of  -5-^. 
_ Sl¬ 
av  8V 
_ x  __  _ s 

By  By 


Differentiating  Eq.  114  we  obtain 

cosaf  -  V  sina  f£.  (123) 

s  ay 


At  the  origin  a  ~  0,  and  this  become. 

av  av 

_ x  =  s 

By  By 


(124) 


But  this  is  merely  the  directional  derivative  of  Vo  in  the  y 
direction,  which  at  the  origin  coincides  with  the  direction 
normal  to  the  line  of  flow.  Equation  124  thus  gives 


8V  8V 
_ x  m  _ s 

8  y  8s„ 
J  n 


(125) 


where  sn  indicates  the  direction  of  the  outward  normal  to  the 
line  of'  flow. 


Deterimination  of 


8V 

x 

unr 


8V 

_ x 

Bz_' 

BVS 

sr 


Differentiating  Eq.  114  we  obtain 

(126) 


coso  -  V 


.  80 

sino 


At  the  origin,  where  or  *  0,  this  becomes 

av  av 

x  ^  _ s 

TT  ~8T! 


(127) 


It  should  be  noted  that  the  direction  of  the  positive  *  axis 
here  is  the  same  as  that  of  the  t  axis  in  Eq.  85.  In  apply¬ 
ing  Eq.  127  we  may  therefore  use  the  s  coordinate  of  Eq.  85. 


ay 

Determination  of  -j~*. 


Differentiating  Eq.  118  we  obtain 


8vy  8V,  »« 

TT  "  '  "!T  sina  '  v.  eo*“  Ti- 


At  the  origin  a  ■  0,  and  this  becomes 


(128) 
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_  v  Ba 
dz  ~  "  s  8z  * 


(129) 


With  Ea.  129,  as  with  Eq.  127,  we  may  use  the  z  coordinate 
of  Eq.  85. 

8V 

Determination  of  j— -.  Substituting  Eq.  117  in  Eq.  113 


we  obtain  immediately 


8Vy 

8y 


9Vs 
8s  * 


(1?0) 


If,  however,  we  differentiate  Eq.  118  we  obtain 
8V  9V 

"5 —  =  -  —a—  sma  -  V  cosa  77-.  (131) 

8  y  <V  s  8y 

At  the  origin  a  =  0,  and  the  y  direction  is  that  of  the  out¬ 
ward  normal,  hence  Eq.  131  becomes 


(132) 


That  Eq.  132,  which  we  have  just  obtained  by  differentiation, 
is  compatibl  'ith  Eq.  130,  which  was  obtained  as  a  conse¬ 
quence  of  the  incompressibility  of  the  fluid,  can  be  seen  by 
examining  an  elementary  length  of  a  tube  of  flow  (see  figure 
below).  Fluid  enters  the  left  face,  of  area  dsndz,  with  a 
velocity  V  ,  and  leaves  the  right  face,  of  area 

w 


(ds 


n 


8a 

8s 

n 


dsnds)dz, 
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with  a  velocity  (Vs+dVs).  Since  fluid  leaves  the  elementary 
volume  at  the  same" rate  at  which  it  enters,  it  follows  that 


V  ds  dz  *  (V  +dV  )  1  -  7p~  ds  ds  dz 

s  n  s  s  \  0sn  J  n 


0  =  dV  -  V 


(133) 


s  0s 


wherein  the  higher  order  term  has  been  omitted.  Equation  133 
may  be  written 


_ s  =  y  Qa 

8s  s  8s  ’ 
n 


(134) 


which  equation  expresses  the  identity  of  Eqs.  130  and  132. 
Equation  130  is  easier  to  apply  than  Eq.  132,  since  its 

8V 

negative  is  already  available  as 

Substituting  Eqs.  117,  122,  125,  127.  129,  and  130  in 
Eq.  43,  we  obtain 


Power  loss  r  /8V  \a  /flV  V  \*  /8V  \a 

ner  unit  =  U  41  ^  -—■)+(  -k-5-  ) 


per  unit  =  M  41 -gj 
volume  *-  v 


©■*(*,«)■]■  <>»> 


We  shall  now  substitute  Eq.  110  in  this  expression,  and 
thereby  determine  the  power  loss  per  unit  volume  good  through 
first-order  terms  in  b. 


Along  a  line  of  flow  A  and  z  are  constant;  hence 


!!s  m 

9s  99  9s’ 


(136) 


But,  noting  Eq.  110, 


[l  ’  ©1 sln9; 


(137) 


also,  noting  Eq.  85  and  applying  the  binomial  theorem, 
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Bfi 

8s 


/a  ,  /'Srla 
^  r  +  (a?) 


=  jra  -  X(ra-a)  -  A*b 


Be 

8s 


/  2  \ 

' "  vv  Jc 
a _ I _ |1+ _ **£-_  [\  fe- 

ra“Mr2-a)  I  ra-X(r2-a)  Vzj 


(138) 


GOS0  + 


•  «  • 


Equations  137  and  138  in  Eq.  136  now  give 


9V 

s 

8s 


2q(A , z )Ab 


z 

z. 


sinfl 


(ra-a)  [r2-A(ra-a) ] 
9V 


Next  let  us  determine 


8s  ’ 
n 


We  have  for  fixed  z 


(139) 


9  V 


9V. 


dvs  “  Tf  dx  +  Tf  d*- 


(140) 


In  the  normal  direction  d0  is  given  in  terms  of  dA  by  Eq.  106 
Substituting  in  Eq.  140  and  dividing  by  -ds  , 4  given  by  Eq. 
107,  we  obtain  n 


8V 

w 


9  A 


8Vs  8r  8r 
Bd  Be  BA 


r  + 


78t\i 
V  B$ 


n 


8r 

9A 


A7!?!}1 


- - - - -  9V 

Here  the  minus  sign  is  due  to  the  fact  that  -5— ^  is  the 

os 

n 

directional  derivative  of  Vg  in  the  outward  direction  normal 
to  a  stream  line . 
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i  +  (±  ®sY 
x  Vr  W 


i  8V  n 

_1 _ s  8r 

a  do  d$ 

jo _ _ 

A  +  Zt 

'  1  +  Vr  80/ 


(141) 


Substituting  Eqs .  85  and  110  in  Eq.  141  we  obtain 

...  ^  ^  -  fe)'M 

n  (r2-a)  jra-a+2Ab|^l  -  (~"\  cosej- 

2q(A,z)  b  1  -  (~\  jcosfl 
(r2-a)a  {ra-a+2Xb  1  -  (j-^j  jcvsflj- 

& a  - . {i  -  ^  i  -  (f)  M 

6sn  (r  -a)*  1  ri'a  L  'V  J  J 

*  (142) 

+  ■iSfei.’Ufe  1  .  (-5-^  Icosfl. 

(r^-a)  *"  Xv2x/ 

Continuing  we  next  use  Eqs,  100  and  110  to  obtain 

Zi  «  - -<T^Hrr - rr  I1  -  [i  -  (7-)  ]cos«j.  (143) 

R  (rj-aHrj-Xfrj-a)]  (  r2-a  |_  Va^/  J  ! 

Next  we  obtain  from  Eq.  110 

ay  r  .  r  ✓  v  n  ^ 


&  - 

8  z  r 


Mr  {*  -  ^  [‘  ■  (t)  >“•} 


(144) 


+  4QXl,z)l_bz  oos(_ 
zi(ri-a? 


Lastly  we  shall  determine  Noting  the  figure  on  the 

following  page  we  see  that 
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a  =  tan 


-1  ft  8r 


8rN\ 

Jej' 


(145) 


8a 

8z 


8ar 

dzdO 


dr  dr 
d0  dz 


I  ,  ( dr 

r  +  U?. 


(146) 


Noting  Eq.  85  we  see  that  the  first  and  second  terms  in  the 
numerator  are  of  first  and  second  order  in  b,  respectively; 
and  that  the  first  and  second  terms  in  the  denominator  are  of 
zeroth  and  second  order  in  b,  respectively.  It  follows  that 


da  a  1  82r 

8z  r  8z00’ 


or,  substituting  from  Eq.  85, 


da 

dz 


3 


2b\az _ 

z* [ra-*(ra-a)  ] 


sin0 . 


This  with  Eq.  110  now  gives  finally 


V 


da 
s  dz 


3 


2Q(X .z)bXaz _ 

zj(ra-a) [ra-X(ra-a)] 


(147) 


(148) 


sin0 . 


(149) 
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We  have  now  determined  all  of  the  quantities  which  appear  in 
Eq.  135,  the  power  loss  per  unit  volume. 

In  order  to  obtain  the  total  power  loss  P  we  must  inte¬ 
grate  Eq.  135  throughout  the  volume  of  the  mercury.  Using 
the  x,  y,  z  axes  of  the  figure  on  page  31  we  have 


x  =  r  cos0,  y  =  r  sin0,  z  =  z  (150) 


where  r  is  given  by  Eq.  85.  The  Jacobian  of  x,  y,  z  with 
respect  to  X,  0,  z  is 


9(*,y,z) 

3(X,0,z) 


f d  r  \  8  r> 

'M  COS9  -  r  sin 6)  ^ 

/  0 


(151) 


Developing  by  minors  with  respect  to  the  elements  of  the 
last  row,  we  obtain 


8.(x,y,z) 

0(X,0,z) 


dr 

dX 

9r 

d\ 


cos  6 


sin0 


cos0  -  r  sin0 


sin0  +  r  cos0 


g.C.X|,y,|Z)  =  £r 
d(X  ,$  yZ)  dX' 


(152) 


The  element  of  volume  in  the  curvilinear  coordinates  X,  0,  z 
is  therefore 


dv  = 


9(x,y,z) 


8  (X  ,  0 ,  z ) 


dXd0dz  ~  *  r  ^  dXd0dz 


(153) 


8r  . 


since  ^  is  negative,  as  is  evident  from  Eq.  85.  Equation  85 
in  Eq.  153  now  gives  finally 

dv  =  |(ra-a)  (rl-X(rJ-a)  1 


Xb  (2ra-3X(ra-a)  j  £l  -  (jt~)  cos 6  (154) 

-  j  cos2©!"  dXd0dz. 


-  2xV 
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Now  that  we  have  obtained  all  of  the  quantities  which 
appear  in  Eq.  135,  and  also  the  element  of  volume,  we  are  in 
a  position  to  write  down  the  expression  for  the  total  power 
loss  P.  Substituting  Eqs.  139,  142,  143,  144,  and  149  in 
Eq.  135,  noting  Eq.  154,  and  integrating  throughout  the 
volume  of  the  mercury,  we  obtain 


pi  p  zi  p  a*  i 

pa"  J.  j 


4qAb 


1  - 


2 
Z  . 


sin0 


(r2-a)  [ra-A(ra-a)] J 


.  8A  f  _  4Ab 
(ra-a)2  ^  r* "a  - 


1  - 


COS0 


I 


2qb 


(ra-a) 


1  - 


COS0 


(155) 


2Ab 
a 


_ 3 _  U  . 

[ica -X Cra-a)  ]  \  ra- 

iS  r 

8z  f,  2A 
ri ”a  I  “  ra“ 


1  -  — 
\z.  j 


COS0 


2Ab 
a 


l  -  f  }  cose r  + 


4qAbzcos0 

2  f  s2 

Zjtrj-a) 


\ 

2<^b\a z  sinfl - 1’(  j(r  ,a)U  -X(r  -a)J 

f r>  -A/'-r  _oMJ  i  l  2  2  2 


Lz*  (ra-a)  [ra-A(ra-a)  ]-*  j 


+  Ab [2ra-3A(ra -a)  ] 


dAdzdtf . 


Here  q  =  q(A,z).  If  we  knew  q  and  b,  Eq.  155  would  give  P 
closely;  however,  q  is  an  unknown  function  which  we  shall  now 
determine.  In  so  doing  Eq.  155  will  be  reduced,  and  the  low¬ 
est  orders  of  terms  in  b  sorted  out.  Further  reduction  at 
this  point,  however,  would  be  premature. 
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DETERMINATION  OF  q(A,z) 


Since  q(A,z)  depends  upon  b,  let  us  expand  this  function 
in  a  Maclaurin  series,  thus 

q(A,z)  =  q0(A,z)  +  bq^AjZ)  +  b2qa(X,z)  +  ...  .  (156) 

The  problem  of  determining  q(A,z)  then  becomes  that  of  obtain* 
ing  q0(A,z),  q1(A,z),  q  (A, z },...,  which  functions  are  inde¬ 
pendent  of  b.  This  we  can  do  by  minimizing  the  viscous 
frictional  power  loss  P,  given  by  Eq.  155,  under  the  condi¬ 
tions  that  the  mercury  velocity  be  that  of  the  channel  walls 
wherever  it  is  in  contact  with  these  walls,  and  that  b  be 
fixed  at  some  arbitrary  small  value. 

First,  let  us  place  b  =  0,  then  Eqs.  156,  110,  and  155 
become,  respectively, 


q(*»z)  =  q0(*>z)> 


Vs  ■ 


q0  (*•*) 


(157) 


(158) 


p„  =  •> 


%  T 


r,-a 


(159) 


(r^-a)  [r|-A(r|-a)  ]dAdzd0. 


Wien  b  =  0  the  mercury  takes  the  shape  of  an  annular  ring, 
and  rotates  with  the  main  spinning  body  as  if  it  were  a 
rigid,  integral  part  of  it,  the  apparent  angular  velocity 
being  (P-v)  clockwise.  Noting  Eq.  85  with  b  *  0,  it  follows 
that 


Vs  *  (P-v)  [rJ-A(rj|-a)  J. 

Substituting  in  Eq.  158  we  now  obtain 

q0(*,z)  “  (P-v)  (ra-a)  fra-A(r|-a)  ]. 


(160) 


(161) 
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Substituting  Eq.  161  in  Eq.  159  we  find  that  both  brackets 
in  the  integrand  of  Eq.  159  vanish,  and 

P0  =  0,  (162) 

which  we  should  expect,  with  the  mercury  rotating  like  a 
rigid  body.  It  may  be  mentioned  that  Eq.  161  can  also  be 
obtained  directly  by  minimizing  Eq.  159  subject  to  the  above- 
mentioned  surface  conditions,  no  use  then  being  made  of  Eqs. 
158  and  160. 

Next,  let  us  substitute  Eq.  156  into  Eq.  155,  which  can 
be  visualized  without  rewriting  Eq.  155;  then  the  resulting 
expression  for  P  becomes  a  power  series  in  b,  thus 

P  *  P  +  P.b  t  P  b2  +  . . .  .  (163) 

0  4  3 

Referring  to  the  subscript  on  or  q^  as  the  order  of  that 

Pi  or  qj_,  je  see  that  the  higher  the  order  of  P-;  ,  the  higher 
will  be  the  order  of  the  highest  order  q  which  it  contains. 

In  order  to  determine  successively  qt(*,z),  qa(*>z),  ...  we 
make  use  of  the  following  two  facts . 

1.  q(X , z ) ,  and  hence  qj(X,2),  qa(X,z),  ...  are  such  that 
P  is  a  minimum,  subject  to  the  above-mentioned  surface  condi¬ 
tion,  regardless  of  the  particular  value  chosen  for  b. 

2.  If  b  is  taken  to  be  sufficiently  small,  the  sum  of 
all  the  terms  ueyond  the  nth  in  Eq.  163  is  negligible  in 
comparison  with  the  sum  of  the  first  n  terms  of  this  series, 
providing  on.t.y  that  these  n  terms  contain  at  least  one  non¬ 
vanishing  tern.. 

We  have  already  seen  how  qf(X,z)  can  be  obtained  by  using 
but  one  term  (n=l)  of  Eq.  163  in  connection  with  the  minimum 
principle . 

In  order  to  determine  let  us  place  n**2;  however, 

since  PQ  contains  only  q0(X,z),  which  has  already  been  deter¬ 
mined,  and  since  PQ  *  0,  due  to  Eq.  162,  we  shall  in  effect 
try  to  determine  qA(X,z)  by  minimizing  Pj ,  in  which  qQ(X,z) 
has  been  replaced  by  its  known  expression  (Eq.  161). 

Substituting  Eq.  156  in  Eq.  155  each  large  bracket  in  the 
integrand  becomes  a  power  series  in  b.  These  series,  however, 
begin  with  terms  which  are  linear  in  b,  for  those  terms  with¬ 
out  b,  which  are  also  those  which  appear  in  the  brackets  in 
Eq.  159,  all  vanish,  as  we  have  seen.  Noting  the  squares  on 
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the  large  brackets  in  the  integrand  of  Eq.  155,  it  is  now 
evident  that  the  smallest  power  of  b  which  can  appear  in  the 
expansion  of  P  is  two;  hence 

?l  =  0.  (164) 


In  view  of  this  we  place  n=3,  and  try  to  determine  qi(X,z) 
by  minimizing  P  . 

Since  the  coefficient  of  b  in  the  series  obtained  by 
squaring  a  power  series  in  b  which  begins  with  the  linear 
term  is  merely  the  square  of  the  coefficient  of  b  in  that 
linear  term,  we  can  write  immediately 


P.  = 


-rpr 

Jo  J-z1  ,;0  I 


4o  A 
‘0 


i  - 

\ZXJ  j 


(ra-a)2  [r1-A(rI-a)  ]  J 


r  9% 

4^X  * 

L(ra-a)S 

9 S 
ax 


-  W1 

v^y  J 


cos$ 


(ra-a)“  (ra-a)  L 


-  (f)} 


cosO 


i  - 


(?' 

KZi 


cosfi 


+  /  \  /  \  i  (T  -a)'fr  -A (r  -a) 1 

(ra-a)  lra-X(ra-a)]  t  *  s  /J 


9,‘o 

2  * 
(ra-a)* 


(165) 


-  (t)  } 


eq 

4q0A*co»* 

cosO  + 


r  -a  i  <  » l 

i  zJ(rJ-a) 


r  2qftA*z  sin®  1*  ) 

— - -0 -  (r  -a)[r  -X(r  -a)jdAdzd*. 

^z*(ra-a) fra-X(ra-a) \ 


Noting  that 
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(A+Bcos0)2dtf  =  jt(2A2+B2)  (166) 


we  square  the  brackets  in  Eq.  165  and  carry  out  the  integra¬ 
tion  on  0;  then  substituting  Eq.  161  for  qQ  we  obtain 


=  tmC/S-p)2 


(r,-a)2 


+  4  [ra-2A(rJ-a)  ]2 


1642z2[f,-X(ra-a)]2  4>y  | 

<  s 


4  /  \  i 

-I  (x*a-a) 


(167) 

(ra-a) fr2-A(ra-a) JdXdz 


+  2*m 


,  Z,  i 

-i 

r 1  r  1 ' 

B\  j,  qi 

(ra-a)2 

*( 

^ (r.-a) (r^-Afr^a)  ]dAdz. 


Here  the  first  integral  does  not  contain  q^  ,  and  hence  plays 
no  role  in  the  minimization  process.  The  integrand  in  the 
second  integral  contains  two  squared  brackets,  and  is  hence 
positive.  This  integral  therefore  takes  its  minimum  value, 
zero,  when 


qj(A,z)  =  0, 


(168) 


which  is  consistent  with  the  required  boundary  conditions, 
and  is  hence  the  desired  function. 


At  first  sight  it  may  seem  strange  that  Pj  *  0  and 
(^(AjZ)  =  0;  however,  noting  Eqs.  85,  110,  and  the  figure  on 
page  31,  we  see  that  if  b  is  replaced  by  (-b),  the  effect  is 
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merely  to  rotate  the  flow  pattern  as  a  whole  through  180  deg, 
which  does  not  alter  P  or  q(A,z).  P  and  q(A,z)  must  there¬ 
fore  be  even  functions  of  b;  hence 

P.  =  0  if  i  is  odd, 

1  (169) 

q^(A,z)  -  0  if  i  is  odd. 

We  note  that  Eqs.  164  and  168  are  in  accord  with  this  result. 

Finally  we  note  that  since  q0(A,z)  provides  the  desired 
mercury  velocity  at  the  channel  walls,  it  follows  that  all 
higher  order  q’s  must  Vcni  i  at  tv  channel  walls,  thus 

q^(A,z)  =  0  when  i  >  0.  and  A  *  0  or  z  =  ±  .  (170) 


DETERMINATION  OF  THE  VISCOUS  FRICTIONAL  POWER  LOSS  P 


Returning  to  Eq.  167  we  see  that  since  the  second  integral 
vanishes  due  to  Eq.  168,  Pa  is  given  by  the  first  alone. 
Carrying  cut  the  integration  on  z,  we  obtain 

j  p  i  r!6Aai6z1  4  (rI-2A(r1-a)  ]ai6zt 

P  =  \  i - r-  + - j - — 

Jo  U5(r  -a)  15(ra-a) 

(171) 

16A*2z  fr  -A(r  -a)]*  4A42z  n 

+ - ---T - +  ~ 7 H  (V*> <Vx<r. ->l* 

3z, (r.-a)  3f,  J 


n  l  r 32A*z 

1Z  .)  Ir«-x(r»-a^ 


8z 


+  - l—  (ra  +  A(-rj(r1-a)-4rj(r1-a)) 


(rJ-a) 


+  A* (4rt (r|-a)*  +  4ri(rI-a)t)  +  A* (-4(ra-a)S ) } 

+  n~^r  u*d-«sd(r.-a)  +  3XVr t -*>* 


l '  i 


.  ,  5(r  -a)  1 

-  xhr,-.)’]  +  - rA*r,-x‘(r,-*i)j- 


dA 
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8  /n  n2  f  32zi  fl  1  , 

=  n  tm(M')  {‘(T^aT  [J  r2  -  4  ra "a 


8z. 


+ 


<Va)s  L  * 


s  5a,  s  ,  8  ,  Na  ,  .3 

r-  '  2  ra'ra"a)  +  3  ra(r2"a)  “  ^ra“a) 


20 


zi(Va) 

£  <va) 


Is  3a,  >.3  ,  v2 

3  r2  -  4  Wa)  +  5  r2(ri-a) 


5(i'a-a)  r 


1 

—  T 

5  ‘a 


6  <ra-a) 


,a  f  ,,  .  32zira  20ziL  d2ira 

r  16z* +  (v?)  -  7^17  -  ^ 

25  (r^-a^4  13r,^r2-a)  15r* 


8z.  r ' 


20  ra 


3  vv*); 


=*}  • 


p  - 


4  5 


-<»->■  {»■.[»  C*)‘  - 5  &)’  •■(*)-  •] 


■"  r»  \ 


7  I-40  W 


90  +  78 


Ta-as 


25 


irfl  ^ 

, )  Jr 


(172) 


Knowing  Pa  and  the  fact  that  P3  =0,  due  to  Ea.  169,  it 
fellows  from  Eq.  163  that  accurate  through  the  term  contain¬ 
ing  b3 


P  ^  P  b2,  (173) 

a 

We  saw  on  pages  15  and  16  that  this  quantity  is  of  primary 
importance  in  eliminating  the  wobble. 


EXTENSION  OF  THE  MINIMUM  PRINCIPLE  FOR  THE  VISCOUS 
FRICTIONAL  POWER  LOSS 


The  class  of  velocity  distributions  considered  in  deriv¬ 
ing  the  principle  of  least  viscous  frictional  power  loss 
(pp.  17-25)  consisted  of  those  distributions  which  at  each 
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point  have  the  right  direction,  and  which  also  satisfy  the 
actual  surface  boundary  condition  at  the  solid  surface,  as 
described  on  page  19. 

Later  on  in  applying  the  minimum  principle  to  determine 
the  velocity  distribution,  the  shape  of  the  velocity  field 
was  specified  by  Eq.  85,  which  was  assumed  to  give  a  suffi¬ 
ciently  good  approximation  if  b  is  properly  chosen.  We  still 
have  the  problem'  of  determining  b. 

It  has  been  convenient  to  picture  the  mercury  flow  as 
being  composed  of  an  indefinitely  large  number  of  elementary 
tubes  of  flow,  a  typical  tube  comprising  the  region  X  to 
A  +  dA  in  A,  and  z  to  z  +  dz  in  z.  We  may,  in  fact,  consider 
the  fluid  to  be  constrained  to  move  in  a  network  of  such 
tubes.  If  in  so  doing  we  picture  the  walls  of  the  tubes  as 
having  the  ability  to  exert  constraining  pressure  normal  to 
these  walls  and  hence  normal  to  the  tubes  of  flow,  we  arrive 
at  a  situation  wherein  the  approximate  velocity  distribution 
obtained  for  any  value  of  b  and  given  by  Eqs.  85,  110,  and 
156  may  be  considered  to  be  the  exact  velocity  distribution 
that  is  obtained  if  a  suitable  constraining  body  force  3  per 
unit  volume  be  superimposed.  The  body  force  intensity  vector 
0  is  everywhere  normal  to  the  relative  velocity  vector  VR. 
Adopting  this  point  of  view  we  shall  henceforth  refer  to  the 
approximate  velocity  distribution  corresponding  to  any  value 
of  b  as  the  "constrained  velocity  distribution . ” 

We  can  easily  extend  the  minimum  principle  given  on  pages 
l7  to  25,  so  that  it  is  applicable  not  only  to  the  actual 
velocity  distribution  but  also  to  the  constrained  velocity 
distribution.  In  order  to  do  this  we  need  only  to  add  a  term 
0  to  the  right-hand  side  of  Eq.  59,  which  leads  to  an  addi¬ 
tional  term 


3*<5Vr  dv  (174) 

on  the  right-hand  side  cf  Eq.  60.  This  term  vanishes,  how¬ 
ever,  since  3  is  normal  to  dVR,  which  is  parallel  to  VR; 
hence  Eq.  60  is  not  altered,  and  the  derivation  continues  as 
before.  It  follows  that  the  principle  of  least  viscous  fric¬ 
tional  power  loss  applies  not  only  to  the  actual  velocity 
distribution,  but  also  to  the  constrained  velocity  distribu¬ 
tion  . 

As  a  consequence  of  this  extension  of  the  minimum  princi¬ 
ple  we  can  dispense  with  the  assumption  that  the  shape  of  the 
velocity  field  given  by  Eq.  95  is  sufficiently  close  to  that 
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of  the  actual  distribution  to  permit  the  application  of  this 
minimum  principle,  since  the  validity  of  this  principle  has 
now  been,  established  for  the  constrained  velocity  distribu¬ 
tion  itselr. 


PRINCIPLE  OF  LEAST  CONSTRAINT 


Since  the  constraining  body  force  intensity  3  was  intro¬ 
duced  artificially  so  that  with  it  the  approximate  velocity 
distribution  becomes  exact,  we  see  that  the  exact  solution 
of  the  fluid  flow  problem  is  characterized  by  the  relation 

3=0.  (175) 

This  will  automatically  be  realized  if 

r* 

\  3*  dv  -  0,  (176) 


the  integration  extending  throughout  the  volume  of  the  mer¬ 
cury.  If  the  solution  obtained  is  not  exact,  the  integral  in 
Eq.  176,  whic’  we  shall  call  the  "constraint,"  will  not  vanish. 
Nevertheless  the  problem  of  obtaining  an  exact  solution  in¬ 
cludes  that  of  obtaining  zero  constraint;  and  in  obtaining 
our  approximate  solution  we  shall  try  to  achieve  this  condi¬ 
tion  as  closely  as  possible,  as  indicated  by  the  condition 


Constraint  =  k 


dv  =  minimum. 


(177) 


At  present  we  have  a  family  of  approximate  solutions 
characterized  by  the  parameter  b  and,  as  we  shall  see,  by  an 
unknown  function  of  A  and  z.  These  will  be  determined  by  the 
above  "principle  of  least  constraint,"  which  is  expressed  uy 
Eq.  177,  and  by  a  surface  analogue  of  this  Drinciple  which 
applies  to  the  surface  constraining  pressure  at  the  iree 
mercury  surface.  This  analogue  will,  however,  be  considered 
later  (see  page  64). 
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DETERMINATION  OF  THE  PRESSURE  FIELD 


If  the  constraining  force  intensity  3  be  included  in 
Eq.  59  we  obtain  as  the  condition  for  dynamic  equilibrium  in 
the  liquid 

PVr*Wr  =  -  Vp  -  pV$  -  2 p/fxVR  +  /iV*VR  +  3.  (178) 

We  shall  use  this  expression  to  determine  3  for  use  in  Eq. 
177;  however,  in  order  to  do  this  we  must  first  determine  p, 
the  pressure. 

p  can  be  partially  determined  by  solving  Eq.  178  for  Vp, 
multiplying  by  lsds,  and  integrating  along  a  line  of  flow 
starting  at  the  point  for  which  0=0,  as  shown  in  the  figure 
below.  Proceeding  thus,  we  obtain 

s  f  s 

ls-Vpds  =  \  [-  Pls'(VR-Wp)  -  pls-V» 

” 0  (179) 

-  .’pls-0XVR  +  |ilsVvR  +  ls-8]ds. 


But  the  triple  product  1s*/SXVr  vanishes  since  1«  and  Vr  have 
the  same  direction;  also  Is*3  vanishes  since  3  is  normal  to 
the  line  of  flow.  Finally,  noting  Eq.  6t  and  the  fact  that 
we  have  denoted  jVR|  by  Vg ,  we  have 
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V(VWR>  =  |VR 


d  |  V 


R 


ds 


1  d<v>,> 

2  ds 


(180) 


In  view  of  the  above,  Eq.  179  becomes 


d(Vl) 


f  8  d£  ds  «  f  L  £L  .  p  d*  +  ul  Vv 

Jn  ds  J  2  ds  p  ds  M±s  VR 

u  n  ^ 


ds , 


PV 


P*- 


-  pi  +  M  j*  ls*VaVRds  +  V>(A,z)  (181) 


wherein  ^  is  constant  for  any  line  of  flow,  but  may  differ 
from  line  to  line.  We  note  that  a  line  of  flow  is  specified 
by  the  pair  of  values  (A,z)  which  correspond  to  it;  hence  $ 
is  a  function  of  A  and  z. 

Taking  the  gradient  of  Eq.  181  gives 

Vp  =  -  pV  W.  -  pV*  +  mV  f  1  •v'Vpds  +  V*.  (182) 

S3  J  q  S  K 

Substituting  in  Eq.  178  and  solving  for  3  we  now  obtain 


3  -  P(Vr’Wr-Vs7Vs)  +  2 P0XVR 

-  »{V*VV  J  *  ls-V*VRds)  +  V*. 


(183) 


Noting  Eq.  68  and  the  fact  that 


8V  8V  8V 

™s‘ljT  h+TTh+T?  V 

n  n 


(184) 


also  that 


1sn  *  -  in-  'V  -  VS! 


(185) 
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it  follows  that 


(186) 


Equation  183  thus  becomes 


/  V*  8V  \  8V 

8  =  <  f-  +  Vs  wj  h  -  pVs  TT  + 

n' 

-  v  j"/  V^V3)  + 


(187) 


Here  we  have  seen  that  3  is  at  each  point  normal  to  the  line 
of  flow.  The  same  is  true  of  V$,  since  ^  is  constant  along 
each  stream  line;  also  it  is  evidently  true  of  the  first 
three  terms  on  the  right-hand  side  of  Eq.  187,  including  the 
Coriolis  term.  It  follows  that  the  one  remaining  term,  which 
contains  n ,  is  at  each  point  normal  to  the  line  of  flow. 

Another  fact  to  be  noted  is  that  since  3  is  by  definition 
such  as  to  make  the  approximate  flow  physically  possible, 
there  is  no  question  about  the  existence  of  hence  the  ex¬ 
pression  for  obtained  by  solving  Eq.  IB''  must  be  such  that 
a  function  having  it  as  a  gradient  exists.  The  curl  of  this 
expression  must  therefore  vanish. 

Turning  to  the  -eterminatic \  of  #  we  first  note  that  Eq, 
187  may  be  written 

3  “  U(0,A,z)  +  V#(A,z)  (188) 

where 


/  v*  av«  \  *V  - 

u  -  <  f +  vs  sr  )  AK  -  *v.  -r#  h  +  *>«VR 

-  <*(?VR  '  V  [S 


(189) 
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The  configuration,  the  velocity  and  pressure  fields,  the 
constraining  force  field,  and  ip  also  depend  upon  the  parame¬ 
ter  b,  which,  in  turn,  depends  upon  {,  which  is  a  measure  of 
the  angle  between  the  spin  axis  and  the  precessior  axis. 
Substituting  Eqs.  188  and  152  in  Eq.  177  we  obtain  for  the 
constraint 


p  l  p  zi  r  a* 

«  J  J  J  (U+V*)  JdAdzdfl  (190) 


-zi  0 


wherein 


J  - 


.  8(x,y,s) 

i8(X,0,z) 


=  -  r 


8r 

81 


{(r^-a)  [rj-XCrj-a)  ] 

Abt^-HC^-a)]  £l  -  f-f) 

1  -  (t)  ]  cos,‘’} 


cos  9 


(191) 


-  21 V 


due  to  Eqs.  153  and  154.  We  now  wish  to  determine  ^  so  as  to 
minimize  x . 

Proceeding  with  ^  and  x  in  a  manner  similar  to  that  used 
with  q  and  P  in  the  section  on  pages  49-53,  we  represent  U, 
ip ,  and  x  by  their  Maclaurin  series  expansions  in  b,  thus 

U(A,z,0)  -  UO(X,Z,0)  +  oUjCX.z,*)  +  b*Ua(X,z,0)  +  ...» 

#(A,z)  «  #0(A,z)  +  +  b%a (X,z)  +  ...» 

,  (192) 

X“X0  +  bXj,tb*#  +  ...  . 

Substituting  in  Eq.  190  we  obtain 

r  ■  x  +bx  +b*x  +...  ■  C  C 
011  Jo  J-z 

i 

+  (Uj-HT^^b  +  (Uj-W^^b*  4-  . . .  ]* JdAdzd# 


60 


NAVWEPS  REPORT  8611 


-  r  i 1  + 2(u«^o) 

0  -Zj  0 

•  )b  +  [(11,+V*,)*  +  2(U0-K^0) 

•  +  ...}  JdXdzd#. 


b  is  evidently  an  odd  function  of  | ,  and  vanishes  when 
1=0,  in  which  case  the  spin  axis  and  the  precession  axis  co¬ 
incide.  For  any  value  of  b,  and  its  corresponding  value  of 
£,  ql»,  and  hence  #0 ,  #a  ,  are  determined  so  that  the 

constraint  k  is  minimized. 

First  let  us  place  b=0,  corresponding  to  which  the  mer¬ 
cury  and  the  main  body  rotate  together  as  one  rigid  body. 
Then  the  second  and  last  terms  in  Eq.  189  vanish,  since  Vs 
does  not  vary  in  the  z  direction,  and  viscosity  plays  no 
role.  Also,  noting  Eqs.  85,  142,  180,  and  161,  we  have 

U0  =  p{(0-*')*  [r,-X(ra-a)  ]  +  (fi-v)*  rra-A(r#-a)]}lN 
-  2p(3(p-v)  [ra-A(ra-a)  ]1^ 

UQ  -  -2p*'(^-v)[r1-A(ra-a)]lN  (194) 

where  In  is  a  unit  vector  directed  toward  the  axis.  Placing 
b«0  in  Eq.  193  we  lose  all  but  the  first  term  in  the  brace; 
and  this  can  be  made  zero  by  putting 

-  V*o  -  U0,  (195) 

thus 

-  v*»  “  3F  ^  lVx<r.-*»V  (1,6> 

d* 

dp 

--JjJ.  *  2py  (0-* )(ra-a)  Ir#-A(ra-a)  J, 

po  -  -  pK/?-*')  (rJ-A(r1-a)]a  4  pv(fi-y)aa  +  c%  (197) 
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where  c«  is  an  integration  constant.  Here  the  second  term 
is  added  artificially  so  that 

*>o  =  c0  when  x  =  1>  (198) 


which  will  be  convenient  later. 

Next,  noting  Eqs.  85,  110,  156,  169,  and  the  figure  on 
page  31,  we  see  that  if  b  is  replaced  by  ^-b)  the  effect  is 
merely  to  rotate  the  flow  pattern  as  a  whole  through  180  deg. 
This  does  not  alter  k ,  which  must  therefore  be  an  even  func¬ 
tion  of  b.  Furthermore  if  we  consider  the  pressures  at 
corresponding  points  of  the  flow  patterns  for  b  and  (-b), 
respectively,  these  being  points  which  differ  only  by  180 
deg  in  0,  the  values  of  A  and  z  being  the  same  for  both;  we 
see  that  all  terms  on  both  sides  of  Eq.  181  except  ip  are  the 
same  at  corresponding  points  (see  NOTE,  below).  It  follows 
that  ip  is  the  same  at  correspondin'*  points.  Since  these 
points  have  the  same  values  of  A  and  z,  but  va  ues  of  b  which 
differ  only  in  sign,  it  follows  that  ip  is  an  even  function  of 
b.  Since  k  and  ip  are  even,  the  k’s  and  tp's  with  odd  sub¬ 
scripts  vanish,  and 

k  =  b  V  +  bV  +  bV  +  ...  ,  (199) 

2  4o 

>(  A,z)  =  ^  o  ( A ,  z. )  +  b%8(A,z)  +  b%4(A,z)  +  ...  .  (200) 


Here  k0=0  due  to  Eq.  193  with  b=0 ,  and  Eq.  195.  It  is  now 
evident  that  good  through  first-order  terms  in  b  the  pressure 
p  is  given  by  replacing  ip  by  Eq.  197  in  Eq.  181,  thus 


pV 


™  -  p  4>  +  p 


y  5  1,-v‘v, 


Is 


-  pi  (0-r)  (r^ -A(ra -a)  ]a  +  py(jJ-n)aa  +  cQ, 


(201) 


NOTE.  In  the  term  p  /0  ls-^VRds  in  Eq.  181  the  inte¬ 
gration  extends  from  the  fixed  point  for  which  0“O  on  a  line 
of  flow  around  clockwise  to  point  s.  The  180-deg  rotation 
of  the  flow  pattern  hence  adds  half  the  line  of  flow  to  the 
path  of  integration  of  the  above  integral  in  determining  its 
value  for  the  point  which  corresponds  to  s  after  the  rota¬ 
tion.  The  integral  will  therefore  be  the  same  for  both  s 
and  this  corresponding  point  if  and  only  if  it  is  not  altered 
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by  this  180-deg  increase  in  its  path  of  integration.  We 
shall  now  show  that  this  condition  is  satisfied. 

Noting  that  «V*VR  is  the  body  force  per  unit  volume  due 
tc  viscous  friction  and  that  this  is  reversed  in  sign  if  VR 
is  reversed  in  direction,  let  us  consider  two  elementary 
lengths  of  the  tube  of  flow  dXdz,  these  being  of  length  ds 
and  located  symmetrically  about  the  line  6=0.  Since  the 
flow  :Ln  one  of  these  is  similar  to  that  which  would  pertain 
in  the  other  if  its  direction  were  reversed  we  see  that 
Mls'^Vj^ds  is  the  same  for  both  of  these  elements  (see  figure 
below).  It  follows  that  the  value  of  nil  -V*  V„ds  taken  over 

S  XV. 


symmetrical 

elements 


the  first  quadrant  0  s  M  f/2  is  the  same  as  its  value  taken 
over  the  fourth  quadrant  3r/ 2  '■  0  X  2t  and  t..at  its  value 
taken  over  the  second  quadrant  is  the  same  as  its  valu*'  taken 
over  the  third  quadrant. 

Turning  to  Eq.  181  we  now  see  th  t  if  we  move  cnce-  around 
the  tube  of  flow  p,  V  ,  and  +  return  to  their  initial  values, 
while  4>  is  constant.  sIt  follows  tnat  the  .ntegrai  must  re¬ 
turn  to  its  initial  value;  hence  the  value  of  Mjij’^VRds 
taken  over  all  four  quadrants  is  zero.  The  above  results  may 
be  written 


W.'W,- 

[ '  i,  *  j,  *  i ,  - » 
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where  the  numbers  indicate  the  quadrants  over  which  the  inte¬ 
grations  are  taken.  Replacing  /  by  /  an^  f  by  /  in  the 
last  e-  "tiori  we  obtain  finally  4  a  3 


Applying  this  result  to  the  term  m  / oS  ls*V*VRds  in  Eq. 

101  we  see  that  since  the  initial  point  of  the  path  of  inte¬ 
gration  is  taken  as  that  for  which  0=0,  it  follows  that  in 
computing  the  value  of  this  integral  for  the  point  which 
corresponds  to  s  after  the  rotation  we  may  omit  the  first 
180  deg  of  the  path  of  integration.  The  integral  over  the 
remainder  of  the  path  duplicates  that  for  point  s,  as  desired. 


SURFACE  ANALOGUE  OF  THE  PRINCIPLE  OF  LEAST 
CONSTRAINT.  DETERMINATION  OF  b  AND  c 

o 


The  constraining  body  force  intensity  2 f  was  introduced 
artificially  to  force  the  mercury  to  flow  in  the  pattern 
specified  by  Eq.  85  (p.  56);  then  the  pressure  was  deter¬ 
mined  so  as  to  minimize  0  in  a  least-square  sense  as  indi¬ 
cated  by  Eq.  177,  which  expresses  the  principle  of  least 
constraint  (pp.  57-64).  In  a  similar  manner  we  shall  now 
introduce  a  "constraining  pressure"  0S  which  acts  on  the  free 
surface  of  the  mercury  from  without,  and  together  with  the 
atmospheric  pressure  pa  and  the  pressure  pt  due  to  surface 
tension  balances  the  pressure  in  the  mercury  and  also  the 
normal  compressive  stress  due  to  viscosity  in  the  mercury  at 
the  free  surface,  as  indicated  by  the  relation 

8V 

3s  -  P  -  2“  -»S  '  pt  -  pa  (202) 


where  n  indicates  the  direction  of  the  outward  (into  the  mer¬ 
cury)  normal  to  the  free  surface.  The  remaining  parameters, 
namely  cQ  and  b  will  now  be  determined  so  as  to  minimize  the 
constraining  pressure  0,.  in  a  least-square  sense  over  the 
free  mercury  surface  a^ ,  as  indicated  by  the  relation 


Surface  constraint 


-£ 


0a  da  =  minimum 


(203) 
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where  da  is  the  element  of  area  of  the  free  mercury  surface. 
This  "principle  of  least  surface  constraint”  is  the  surface 
analogue  of  the  principle  of  least  constraint  (Eq.  177), 
which  was  used  in  the  previous  section. 

Before  proceeding  with  the  determination  of  c0  and  b  it 
may  be  well  to  consider  in  retrospect  the  different  varia¬ 
tional  processes  that  we  have  used.  These  may  be  outlined 
as  follows: 


1.  For  arbitrarily  chosen  values  of  £  and  b,  the  shape 
of  the  relative  velocity  field  is  specified  by  Eq.  85,  after 
which  the  determination  of  this  field  is  completed  using  the 
principle  of  least  viscous  frictional  power  loss. 

2.  The  pressure  field  is  next  determined  except  for  the 
additive  constant  co  using  the  nrinciple  of  least  constraint. 

3.  Holding  |  fixed  \w  finally  determine  cQ  and  b  by  the 
principle  of  least  surf ace  constraint. 

In  the  previous  scccion  the  value  of  b  used  in  determining 
the  pressure  field  was  the  actual  value  determined  by  3 
above,  althouch  at  the  time  we  did  not  know  b  explicitly  in 
terms  of  £ . 


APPLICATION  OF  THE  PRINCIPLE  OF  LEAST  SURFACE 
CONSTRAINT  TO  THE  DETERMINATION  OF  b  AND  c 

o 

In  order  to  carry  out  the  procedure  described  in  tbe  pre¬ 
vious  section  we  must  obtain  expressions  for  the  various 
terms  which  compose  3S  and  hence  p,  which  are  given  by  Bqs . 
202  and  201,  respectively.  We  shall  now  consider  these  terms 
in  order  as  follows: 

pVs 

The  Term  y~.  This  occurs  in  p,  and  follows  from  Eqs. 
110,  156,  and  169 ,  thus 

r  _4xb 

2  2<i'J-a)*  1  va 

which  expression  is  good  through  first-order  terms  in  b. 


-  ft)  r4 


(204) 


65 


NAVWEPS  REPORT  8611 


The  Term  p$„  This  appears  in  p,  and  is  given  by  Eq.  7, 
thus,  noting  the  figure  above, 


p*»  -§  rV  =  -  -P-|-  [ra  +  4*  +  2rs£  cos(fl-0)].  (205) 


Here  r  is  given  by  Eq.  85  with  X^i,  thus 


r  ~  a 
s 


•  ■  G^)1x's9 


(206) 


This  in  Eq.  205  gives,  accurate  through  first-order  terms  in 

b, 


**  “  -  PJ2 -  K  ■  2ab[~  -  (fj)  ]cosa  +  ;S 

+  2|  -ja  -  b£l  -  cosej  cos(0-<£)}, 


p4>  -  -  |a*  +  |a  +  2|a  cos (Q~<p) 


+  b  cosfl 


.L 


-  ©  ]  / 


2a  -  24  cocw-y) 


■, n(i7) 


j- 
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Omitting  the  higher  order  term  which  contains  this  becomes 


=  -  e4a-  +  pa;?1 


■  ©■ 


COS0 


-  k  COS 


mV 


(208) 


Here  it  is  likely  that  the  last  term  in  Eq.  208  is  of  higher 
order  and  can  he  omitted;  but  this  can  be  decided  later. 

f  s 

The  Term  p  \  1  *V“Vr,ds.  This  occurs  in  p,  and  can  be 

_ Jo  s  K 

conveniently  determined  using  the  notation  and  cnoice  of 
axes  described  in  the  section  on  pages  38-48  and  shown  in 
the  figures  on  pages  38  and  39.  Proceeding  thus  we  obtain 


VR  -  ivx  +  jv  +  kvz,  18 


v’v-  =  iVJV  +  jV*V  +  kV*V 

R  X  J  y  5 


is  VvR  -  v\  - 


a2v  a2v 


dx 


r  + 


8y 


i  at  the  origin, 


a2  V 


(209) 


and  must  now  determine  the  three  terms  which  compose  this 
Laplaeean.  The  values  of  these  three  terms  are,  however, 
required  only  at  the  origin. 


Differentiating  Eq.  114  now  gives 


8V 
_ x 

3x 


cosa  - 


sina 


da 

8x> 


02V 

s 

- -  cosa 

8x 


8V 
_ s 

8x 


sina 


8a 

8x 


sina 


8a 

8x 


-  V  cosa 
s 


sina 


il® 

2* 

dx 


(210) 


At  the  origin,  which  is  the  only  quint  in  which  we  are 
interested,  o^O,  and  Eq.  210  becomes 
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82  V  92  V 

_ x  _  _ s 

9x*  9x* 


Similarly,  at  the  origin 

82  V  82  V 
_ x  s 

9  y 2  9y2 


V 

s 


92 V  92  V 
_ x  _  s 

2  2 

9z  9z 


[dzj  ‘ 


(211) 


(212) 


(213) 


In  Eqs.  211,  212,  and  213  Vs  is  given  by  Eq.  11C;  and  the 
following  directional  derivatives 

9a  _  da  da  _  da  9a  .  „ 

9x  8s’  9y  8s  *  8z’  K 

J  n 


which  are  evaluated  at  the  origin,  are  available  from  Eqs 


92V£ 


is  available 


121,  134,  and  148,  respectively.  Also 

from  Eq.  110.  Our  problem  hence  consists  mainly  of  deter- 
9aV  92 V 

'  s 

mining 


8x 


?  and  T 


S 

TT- 


At  the  origin 


(215) 


since  both  sides  of  this  equation  express  the  directional 
8V 

derivative  of  — g"  in  the  direction  of  a  line  of  flow.  Noting 

that  the  direction  of  increasing  sR  is  outward  and  normal  to 
the  line  of  flow,  we  now  have 


8V 

s 

8x 


8V 


s  8s 


8V  9s 
s  n 


9s  8x  9s_  8x’ 


n 


(216) 
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which,  noting  figure  above,  may  be  written 


8V  8V  8V 

=  Ts  COSO!  +  8s^ 


sina. 


(  217  ) 


It  follows  that 


±/dV 

8s  V  0x, 


8aV 


r  cosa  - 


8s 

8a  y 


8Vs  ,  8« 

8S  Sin“  8s 


8V. 


a  .  .  s  8a 

8W  Sln“  +  8i"  cosa  81- 

n  n 


At  the  origin,  where  a=0,  this  becomes,  noting  Eq.  215, 

8V. 


83V  9*  V  ~v  - 

_ s  _ _ s  .  _ s  8a 

a  *  ft  8  9s  8s* 

9x  9s  n 


(218) 


In  a  similar  manner  we  have  at  the  origin 

i .  •  r 

8yJ  Ssn  1 


8*  V 


8y/’ 


(219) 
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since  both  sides  of  this  equation  express  the  directional 
8Vs 

derivative  of  —5—  in  the  direction  of  the  outward  normal  to 

8y 

the  line  of  flow. 


9V 


s  9s 
9s  9y 


8sn 
8y  * 


(220) 


or,  noting  the  figure  above, 


9V 
_ s 

9y 


sina  + 


cosa, 


(221) 


At  the  origin, 


where  o=0,  this 

8aV  9a V 
_ s  _  s 

8y3  9s* 

*  r\ 


become^,  nocing  Eq. 


8V 
_ s 

"8s 


9  a 
8sn* 


219, 


(222) 
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Equations  218  and  222  in  Eqs.  211  and  212,  respectively,  now 
give 


!\ .  !\  +  !!s  u,  „  /a«v 

8xa  8sa  8sn  8s  s  ’ 


(223) 


e\  _  !!s_8«  v  /W\a 

8ya  8sJ  8s  8sn  s  W  ' 


(224) 


Continuing  with  Eq.  223  we  note  that  — g-j  is  given  as  a 

function  of  X,  z,  and  0  by  Eq.  139.  Along  a  line  of  flow  X 
and  z  are  constant;  hence 


_s  =  JL  (  s)  80 
a  8$  V  8s )  8s* 


Applying  Eqs.  138  and  139  this  becomes,  good  through  first- 
order  terms  in  b, 


,  2q(X,z)Xbj^l  -  Jcos0 

S  eg  ^ 

■>J  (ra-a)a  [ra-X(r|-a)  ]* 


(225) 


or,  replacing  q(X,z)  by  Eq.  156  and  hence,  in  effect  by  Eq. 
161, 


8»v  2Xb(/3-^ )  1  -  [f-J  cos 0 

9S*  “  (fj-ajfrj-ACrj-a)]  * 


(226) 


The  second  term  of  Eq.  223  is  given  by  Eq.  142  in  which 
q(X,z)  is  replaced  by  Eq.  156  and  hence  Eq.  161,  and  by  Eq. 

121,  which  is  and  Eq.  100;  thus 
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av 


s  da 


Qs  Bs 
n 


tf-.o.  X 

,-a)  l 


r -Mr 

a  K  2 


1  - 


4Ab 

ra"a 


i  -  U) ) 


cos  0 


+  2(M)b 
(r2-a)2 


1  -  TT- 


(227) 


COS0  . 


The  third  term  of  Eq.  223  follows  from  Eq.  110,  in  which 
q(A,z)  is  replaced  by  Eq.  156  and  hence  Eq.  161,  and  from 

Eq.  121,  which  is  and  Eq.  100;  thus 

V  ($2)*  s  - LP-V±—  U  _  L  _  ^V\  ]C0SA  (228) 

s  \dxj  ra-A(ra-a)  l  ra-a  \z  J  J 


Substituting  Eqs.  226,  227,  and  228  in  Eq.  223,  we  obtain 
finally 


aav 


x  *  2(ff-y )b 


8x3  (ra-a)a 


j.  - 


COS0 


(229) 


good  cl. rough  first-order  terms  in  b. 


av 


N^xt,  continuing  with  Eq.  224  we  note  that  is  given 

8aV  n 

by  Eq.  142;  also  - 1  can  be  obtained  by  replacing  V  by 

8s  S 

0V  n 

in  Eqs.  140  through  141.  Since,  due  to  Eq.  85, 
n 


1  + 


fld£ 

Vr  30, 


1. 


we  thus  obtain 


8aV 


i.  £Z&\ 


0A 


x8s  j 

K  TV 


a  s 


n 


H 

ax 


a  I  80  \3s  J 

r>  L  rK 


n'  -1 


8r 

80* 


(230) 
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Substituting  Eqs.  85  and  142  in  Eq.  230  we  see  that  the 
second  term  in  Eq.  230  is  of  second  order  in  b,  and  can 
hence  be  omitted.  In  the  first  term  q(A,z)  is  replaced  by 
Eq.  156  and  hence  Eq.  161,  and  we  obtain  finally 


8av 


BiP-v  )b 
(r2-a)a 


1 


COS0  . 


(231) 


Noting  Eq.  214  we  substitute  Eq.  231  and  the  expression 
9  c 

for  g-g--  given  by  Eq.  134  in  Eq.  224,  and  thus  obtain 


a 


+ 


6(£-Qb 

(ra-a)a 


(232) 


From  Eqs.  110  and  139,  however,  we  see  that  the  first  term 
on  the  right-hand  side  of  Eq.  232  is  of  second  order  in  b; 
hence  good  through  first-order  terms  in  b  we  have 


9aV 


8y 


X  *  t  , 

a  t  .a  I  x 
(ra-a)  L 


r>. 


cos  0 


(233) 


Turning  finally  to  Eq.  213  we  replace  q(A,z)  in  Eq.  110 
by  Eo.  156,  and  hence  Eq.  161.  Differentiating  twice  with 
respect  to  z  and  placing  A~l,  corresponding  to  the  inner 
free  surface,  we  then  obtain  for  the  first  term  in  Eq.  213 

— s*  cos 0.  (234) 

Bz  (r|-a)zJ 


Next  we  see  from  Eqs.  110  and  148  that  the  second  term  in 
Eq.  213  is  of  second  order  in  fc,  and  can  hence  be  omitted. 
Good  through  first-order  terms  in  b  we  thus  have 

— 3.  a  4.ab.(£t±I  COS0.  (235) 

8z*  (r|-a)z1 
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Putting  A=i  in  Eq.  85  we  obtain 


ds 


r  66 


ds  s 


(236) 


for  the  element  of  a  flow  line  in  the  free  surface.  Substi¬ 
tuting  Eqs.  229,  233,  and  235  in  the  last  equation  of  Eq.  209; 
multiplying  by  Eq.  236  and  integrating,  we  obtain  finally 


M 


1  •VaVDds  as 
s  R 


to’btfep  sin9 

(ra-a)z1 


Qnab(p-v) 

(ra-a)2 


(237) 


good  through  first-order  terms  in  b. 

»V 

The  Term  2m  ”9^*  This  occurs  in  At  any  point  of  the 

free  mercury  surface  let  ln  be  a  fixed  unit  normal  vector 
pointing  outward  (into  the  mercury).  Then 


V  1 


s  s 


1 

n 


since  VD  *  V  1,.  It  follows  that 

K  S  S 


8V  0V 
n  _  s  , 

TH  "  'TFn  Xs 


(238) 


(239) 


It  is  necessary  to  determine  this  quantity  only  on  the  free 
surface.  On  this  surface  ln  t  1^;  hence  the  first  term 
vanishes.  In  the  second  term  we'note  that  along  the  normal 
ln  remains  constant,  whereas  in  general  1  does  not.  Eq.  239 
thus  b  comes 
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av  ai 

_ n  —  w  i  .  _ jl 

9n  s  n  an' 


(240) 


In  order  to  obtain  this  dot  product  we  shall  determine  compo¬ 
nents  along  the  directions  of  the  mutually  orthogonal  vectors 

lc  ,  1-  ,  1„ ,  these  vectors  being  evaluated  on  the  free  sur- 
^o  bn  * 

face  at  the  base  of  the  normal.  Here  the  subscript  o  has 
been  added  to  s  to  avoid  confusing  the  fixed  vector  ls  with 

the  variable  vector  lc ,  which  at  the  free  surface  becomes  lr  , 
as  shown,  below.  ^  ~o 


Although  ls  is  variable,  its  component  along  lz  is  everywhere 
zero;  furthermore  it  is  evident  that  dl,ils  .  It  follows 

81  « 

that  the  only  nonvanishing  component  of  is  that  along  . 

In  Eq.  240  we  therefore  need  'nly  that  component  of  ln  which 
lies  along  1^.  The  component  along  1Sq  is,  of  course,  zero, 

since  1„  is  tangent  to  the  free  surface. 


Noting  the  following  figure  and  the  fact  that  1  41  w* 
see  that  o 


Component  of  ,  CQS 
ln  along  ls 


(241) 
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channel 


mercury 


Moving  along  dsn  from  the  free  surface  to  the  lower  end  of  dz 
in  the  figure  above,  r  and  0  change  by  amounts  dr  and  d0  given 
in  terms  of  dX  by  Eq.  106.  Starting  at  the  same  point  and 
moving  along  the  free  surface  to  the  upper  end  of  dz  we  have 
A=l,  dX-0 ,  and 


dr 


8r 


dz  + 


(242) 


r  being  given  by  Eq.  85.  Since  dz  is  vertical  dr  and  d0  are 
the  same  as  before;  hence  Eq.  242  serves  to  determine  dz. 
Substir  rting  Eq.  106  in  Eq.  242,  we  obtain 


2  3r 
r  SX 


dz  + 


9r 

80 


8r  8r 
80  8 A 


4 


dX 


8r 

8X 


dX 


8r 

8z 


dz 


dz  ~ 


Or 

8X 

8r 

9z 


dX 


(243) 


where  dX  is  the  change  in  X  corresponding  to  dsn.  Equations 
107  and  243  in  Eq.  241  now  give 


Component  of 

1  along  lr 
n 


1 


8r\a 

ddj 


+ 
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Component  of 

1  along  1 
n  3  s 


(244) 


Noting  Eq.  85  it  now  follows  that 


Component  of  1  along  1  si 

n 


(245) 


good  through  first-order  terms  in  b. 

Next  we  shall  determine  the  component  of  1  along  1 ,  . 
Noting  the  figure  on  page  75  we  see  that  s  Sn 

Component  of  1  along  1  =  -  sina;  (246) 

s  ?n 


however,  since  or=0  at  the  free  surface  we  have  in  effect 


a 


da? 


da 

9sn 


(ds  )  i 

n  normal 


da 

dz 


(dz) 


normal 


(247) 


where  the  word  ’’normal”  has  been  added  to  stress  the  fact 
that  these  quantities  correspond  to  an  element  dn  of  the 
normal  line,  and  not  to  an  element  of  length  in  the  free 
surface,  as  shown  in  the  figure  on  page  76.  However,  if  the 
triangle  of  legs  dsn  and  ds  in  the  figure  on  page  76  be 
rotated  90  deg  clockwise,  we  see  that  the  quantities  dz  and 
(-ds«)  in  the  figure  on  page  76  can  play  the  role  of 
(dsn)normal  and  (dz 'normal*  respectively,  in  Eq.  247,  in 
which  case,  noting  Eqs.  107  and  243, 


dr 

d\ 


dX 


dn 


In 

.ex 

8r 

dz 


dX  a 


dX, 


(248) 
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or,  noting  Eq,  243, 

dn  ^  Ids! .  (249) 

Equation  247  now  gives 

da  da  I  dz  I  9_a  ^sn 

dn  9s  dn  ’  9z  dn*  ^  ' 

n 

wherein  the  +  sign  would  be  obtained  if  ln  sloped  upward  in 
the  figure  on  page  76,  in  which  case  the  triangle  would  be 
rotated  90  deg  counterclockwise.  Noting  Eqs.  107,  248,  and 
249  this  becomes 


da  ^  9a  9 a  9r 

dn  ~  9s  9z  9z‘ 

n 


(251) 


From  Eqs.  85  and  149  we  now  see  that  the  second  term  is  of 
second  order  in  b >  and  can  hence  be  omitted.  The  first  term 
is  given  by  Eq.  134,  thus 


da  ~  _1_ 
dn  ”  V 

s 


(252) 


Since  from  Eqs.  246  and  247 


Component  of 


dl 

s 

dn 


along 


sina  _  sinda  _  _  da 
dn  dn  "  dn 


(253) 


at  the  free  surface,  it  now  follows  from  Eqs.  240,  245,  246, 
252,  and  253  that 


9V  9V 

_ n  _ s 

9n  ~  "  9s # 


(254) 


Substituting  Eq,  139  in  which  q(A,z)  has  been  replaced  by 
Eq,  J56  and  hence  Eq.  161,  we  obtain  finally 


2 M 


_ n  fyMP.’.yl 

9n  r^-a 


sin0. 


(255) 


78 


NA7WEPS  REPORT  8611 


The  Term  p  .  This  occurs  in  0  .  Using  the  lines  of 

L  O 

curvature  as  parametric  lines  let  us  consider  an  element  of 
area  of  the  free  mercury  surface,  shown  in  the  figure  below. 


Denoting  the  surface  tension  by  T,  this  being  the  force  per 
unit  length  transmitted  across  any  line  in  the  surface,  the 
normal  component  of  the  force  exerted  by  surface  tension 
along  the  sides  of  length  dsi  is 

Normal  component  ds 

of  forces  on  sides  ~  T  ds  (256) 

of  length  ds^  a 

where  is  a  principal  radius  of  curvature.  Here  the 
reasoning  is  essentially  the  same  as  in  the  case  of  a  vibra¬ 
ting  string  or  membrane.  Adding  to  Eq.  256  a  similar  expres¬ 
sion  for  the  normal  component  of  the  forces  transmitted  along 
the  sides  of  length  dsg ,  then  dividing  by  the  area  dSjdSj  we 
obtain  for  the  pressure  due  to  surface  tension 


6  The  effect  of  surface  tension  is  considered  here  in 
connection  with  the  free  surface  given  by  Eq.  85  with  A*l; 
however,  surface  tension  was  not  used  directly  in  determin¬ 
ing  this  surface. 
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Pt 


=  T 


2T® 


(257) 


where  ®  is  the  mean  curvature  of  the  surface. 


In  order  to  compute  the  mean  curvature  of  the  free  mer¬ 
cury  surface  we  place  A=1  in  Eq.  85,  use  z  and  0  as  parameters, 
and  have  as  the  equations  of  the  free  surface 


x  =  r  cos0  =  a  cos0 


y  =  r  sin0  =  a  sin0 


z  =  z. 


(258) 


Differentiating  we  obtain 


9x  _  2bz 
0z 


2 


COS  0  , 


ay  2bz  .  .  . 

S7  =  — 7  cos0s m0 , 


3z 

8z 

8z 

8x 

00 


=  i; 


=  -  a  sin0  +  2b 


(259) 


|—  *  a  cos 0  +  b 
=  0 

80  u* 


[‘  •  © 
[‘  ■  © 


cos0sin0, 


(3ina0-cos20) , 


Accurate  through  first-order  verms  in  b  the  quantities  E,  F, 
and  G  in  the  first  funoamental  form  are  hence1 


*Struik,  Dirk  Jan.  Lectures  on  Classical  Differential 
Geometry.  Cambridge,  Mass.,  Addison -Wesley  Press,  1950,  pp.  58, 
75,  and  83. 
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Esi, 


F  s  (-cos20sin0  +  cos20sin0)  =  0, 

zi 


(260) 


G  s  a  +  2ab 


-  (-7-)  I ( - 2  sin20cos0  +  sin20cos0  -  co s**0) 

V  Z  j 


G  s  a 


r 

-  2ab  1  -  Jcos0; 


also 


2  2  !  f  7  ' 

EG  -  F  «  G  *  a  -  2abj^l  - 


cos0,  (261) 


from  which  we  obtain  as  the  element  of  surface  area 


daf  =y  EG-F2  dzd0  a  |a  -  b  1  -  (~-^) 


cos0j-dzd0.  (262) 


Taking  second  derivatives  we  now  obtain 

a3x  2b  % 

A  a  2  cos  6> 

9  z  z. 


~ ^  ”  cos0sin0, 

9za  z\ 


9az 

9z2 

9ax 

9z90 


=  0; 


4bz 

a 

*T 

1 


cos0sin0 , 


9ay  2bz  /  a  *  .  a„N 

me =  — r (cos  «-sln  e)' 


(263) 
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a3z 

8z80 


=  0; 


82 

— 7  =  -  a  cos0  +  2b] 
902 


1  - 


z 

z. 


(cos2  0-sin2  0) , 


9a  y 

— 7  =  -  a  sin0  +  4b 
80 2 


1  -  — 


z 

z 


sin0cos0 , 


32  z 

_  -  o. 

80 


Using  Eqs.  259  and  263  we  now  set  up  the  third-order  deter¬ 
minants  which  give  e,  f,  and  g  of  the  second  fundamental  form. 
In  each  case  the  last  column  is  0,  1,  0;  hence  developing  by 
minors  with  respect  to  this  column,  and  omitting  terms  which 
contribute  only  terms  of  order  two  or  higher  in  b  to  the  ex¬ 
pansion  of  tht  determinant,  we  obtain 


e  s  - 


2b  a* 
—  cos  0 


-  a  sin0 


2b 


cos0sin0 


a  cos0 


+  X/EG-F2 


cos0  +  v/EG-F* , 

zi 


(264) 


f  a  - 


~  cos0sin0  Qj  (cos20-sin20) 
z!  zi 


a  sin0 


a  cos0 


f  a  2a-^£  sin0  (2cos20-cos20+sin20) 
zi 


+  /EG-F* 


f  a  sin0  +  y/ZG-?* , 

zi 


(265) 


82 


-  a  cos 6 


+  2bL1  ’  ©*]( 


2  2 

cos  0-sin  0)  -  a  sin0  +  4b 


-  a  sin0  +  2b 


\x  .  (*-)'] 
L  J 


cos0sin0 


a  cos#  + 


+  ""I1 '  ©} 
^  ’  fe)  ](si’ 


+  v/EG-F1 


g  s  -|aa  -  6abj^l  -  sina0cos 0  -  3ab  1  -  (cos* 0- sin* 0) cos flj-  +  >/eg 

g  =  jaa  -  3ab  1  -  (jr^j  cos0(2sina0+cosa  0-sina3)j-  +  v/EG-F^ 
g  ^  -|aa  -  3ab|\  -  cosoj-  *  v/EG-F*  . 


The  mean  curvature  is  now  given  by 


gE-2fF+eG 
2(EG-F*)  » 


aa-3ab  1  -  ]cos0  -  cos0 

2  ja*-2ab  1  -  cos*}  ^ 

i  (0  *  1  (0}°"}  - 

4  f*  -  *S  — )  • 


' 

V  J 


zi 


Substituting  in  Eq.  257  we  now  obtain  finally 


Pt  «  -  ^  (l  -  cos  0^ 


A 
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2b 


2b 


1  - 


1  - 


/V\  , 

\ZJ  J 

f , 


(cos*  0 -sin2 


cos0sin0 


-  a  sin9  +  4b[l  -  j 


sin0cos0 


a  cos0  +  b 


1  - 

L  W 


(sin*0-cos3  0) 


sin  6cos8  -  3ab 


al 


1  “  ( ~z~ )  j(cos30-sina0)cos0j  +  v/eG-F* 


.21 


cos0(2sin20+cos20-sin20)i  +  y/EG-F* 


jz 

2. 


cos0j“  +  /EG-F*  . 


are  is  now  given  by 


gE-  2f  F+eG 
2(EG-Fa )  » 


a  -3ab 


1  - 


2a3b 

COS0  -  -■  T  COS0 


{a*- 


2ab 


cos0r 


1  a/ 2 


J 


COS0 


h  *  ft  ■  ©']-} 


257  we  now  obtain  finally 


pt  “  "  a  (X  ’  oosa)  * 


(266) 


(267) 


(268) 


(269) 


6 
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Placing  Eqs.  204,  208,  and  237  in  Eq.  201;  and  then  plac¬ 
ing  Eqs.  201,  255,  and  269  in  Eq.  202,  we  now  obtain 


>,  *  -  {>  •  ^  [•  -  (f/M 


*  t£A  -  ..»■  {: 


r 


\  s  i 


par  ibl  1  -  {j-)  icos® 


-  (  cos(9-*)}  +  -  S&&1&±  slne 

J  (r2-a)  zx 


_  8/xab(ft-tQ 
(r  -a)2 


.(•X 

Iz,  / 


+  c  .  tokikx  I 

0  ra-a 


sin® 


ai 


uJ 


1  -  (~J  I  sin® 


+  !  (i  - 

Z1 


A. 


(270) 


accurate  through  first-order  terms  in  b.  This  may  be  written 


2L  =  A  +  B  nin®  +  C  cos® 
s 


(271) 


where 


a  a 


a  =  -  tso&ni  +  +  eiil  +  c  +  2  .  p 

2  2  2  o  a  *a 

A  *  %  (2a*^-a*va+^15a)  +  ~  -  p,  +  cn  , 


(272) 


B  «  pa/3  i  sin <t> 


/  \*“1 
f  z  \  J 

\25.y  * 


-  “Si  >'.*•>  [‘  *  (t)’ 


(273) 
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C  —  p  5  b 

i 1  -  fr)  1 

f  Mkall  .  „*1 

i 

L  **  j  ^ 

+  pa0a|  cos*  - 

2Tb 

2  * 

(274) 


If,  noting  Eq.  262,  we  also  place 


D  « 


(275) 


then  Eq. 


K 


S 


203  becomes 

(A  +  B  sin0  +  C  cosiJ)1 

minimum 


(a  -  D  cos0)dzd0 

(276) 


where  ks  is  the  surface  constraint.  Multiplying  out  the 
first  factor  and  carrying  out  the  integration  on  0  now  gives 


(a(2A,+B,+Ca)  -  2  A  C  Djdz 


minimum.  (277) 


Next  we  note  that  in  Eq.  39  P  is  given  by  Eq.  173  accurate 
through  third-order  terms  in  b;  and,  noting  Eq.  92  and  the 
figure  on  page  66,  (-y0)  is  Eq.  93  times  sin*,  accurate 
through  second-order  terms  in  b.  Equation  39  thus  becomes 


P  b*  .  auVtf-M  Sin* 

*  J(rj-a* ) 


wherein  both  sides  are  accurate  through  second-order  terms  in 
b.  Later  we  shall  see  that  due  to  Eq,  291  they  are  accurate 
through  third-order  terms  in  b.  Substituting  from  Eqs.  83 
and  172  tnis  becomes 


/?*4  sin*  *  K  M(/3-^)b 


(279) 
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where 


From  Eq.  279  it  follows  that 


cos <f>  =  ±  7 /34 Z 1  •  (281) 

toe  regard  |  as  given,  whereupon  0  is  determined  by  Eq. 
279,  which  can  be  used  to  eliminate  <jfr  from  Eqs.  273  and  274. 
cQ  and  b  are  now  determined  so  as  to  minimize  the  surface 
constraint  *s,  given  by  Eq.  277.  Let  us  first  determine  c0. 
This  we  can  do  by  minimizing  *  for  an  arbitrary  choice  of  b, 

f-d 

which  amounts  to  equaling  to  zero,  *s  being  regarded  as 

o 

a  function  of  c  and  b.  Noting  Eqs.  272  through  275  we  thus 
obtain  0 

l  _ 

( 2aA- C  D)dz  *  0  (282) 

l 


as  the  relation  which  determines  c0  as  a  function  of  b. 
Multiplying  Eq.  282  by  2 which  does  not  contain  z,  and 
subtracting  fiom  Lq.  277  we  then  obtain 


(-2A  a+B*-£a)dz 


minimum 


(283) 


as  the  relation  for  determining  b. 

Next,  noting  Eqs.  272  through  275  we  carry  out  the  inte¬ 
gration  in  Eq.  2C2,  and  thereby  obtain 
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16 


a 


4az  A  -  z  pab3 

i  15  1  I  r  -u 

*-  a 


y  z^ab/?*^  cos*  +  ~  zj  •— ~  *  0 


A  -  4 


±  pb»  -  *a  |  4  I  PW*{  cos*  -  .  (284) 

L  *  J  3az_ 


This  expression  equated  to  Eq.  272  determines  c  ,  which 
appears  in  J.  alone.  Substituting  Eq.  261  in  Eq?  284  we  now 
get 


*-»*&  4^  •'■]•!  A) 

•  dZ . 


(285) 


±  y  pb 


cvTT^Tv  . 


Equations  279  and  281  in  Eqs .  273  and  274  now  give,  respec¬ 
tively, 


B  *  bMf/?- 


-  •££[■♦  (0  Ctf  -  (0] } . 


C  «  b 


H*  -  (0 1  ■  »•]  •  fi 


*  pa 


4.1  v*  «  ,  ,  .1.1 

(  «K  |i  (P->>)  b  . 


(287) 


Differentiating  with  respect  to  b  the  condition  ",\t  283 
becomes 


(288) 
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Here  A,  B,  and  C  are  given  by  Eqs.  285,  286,  and  287,  respec¬ 
tively;  and  t,  is  constant.  Equation  288  determines  b  as  a 
function  of  £.  Noting  that  does  not  contain  2.  Eq.  288  be¬ 
comes 


-  4z 


b» p  .  A  .1 

i|_  1.15  L  ra“a  J  3 


±  a  Pb  vVi’-kVw-joV 


“fe  '[^ 

"  L  a  J  az. 


2  _T_ 

3  a 

az  - 

l 


/ 

*  i  0  )  - tVfiiifcl! - ) 


( 


vT?-KV{M)‘b 


a.  a 


,  a  p  Z1  r  4(r  +a)  r  / , \>  /r  -a\ 

*  *  i-,.  {'“  •  ^  [*  *  y  y=) 


(289) 


z_ 

z. 


f«  *  C  [4*1'  •  (0} 


2atf-y)  „  „*]  .  2T\ 

J 

J  V 


ra-a 


pa  //jV-kVc^-ioV  pa  1  -  -  2*] 

2T  y  - j  d2  ,  o. 

zi  y$Atx °  J 


X"|*  (fl-v)  b 


This  expression  determines  b  as  a  function  of  £.  Lev  the 
Maclaurirt  series  expansion  of  thrs  function  be 


b  “  b^  +  bsl*  +  b^f*  i-  ... 


( 290 ) 


wherein  the  coefficients  blt  bJ}  ...  are  unknown  constants, 
and  the  terms  containing  even  powers  of  l  are  omitted  because 
b  is  evidently  an  odd  function  of  | .  Good  through  second- 
order  terms  in  i  Eq.  290  becomes 


88 


NAVWEP3  REPORT  8611 


b  ^  . 


(291) 


Substituting  in  Eq.  289  we  see  that  the  first  term  is  of 
third  order  in  whereas  the  other  two  terms  are  of  first 
order.  We  therefore  omit  the  first  term,  the  resulting  ex¬ 
pression  being  good  through  second-order  terms  in  £.  Divid¬ 
ing  this  by  |  we  obtain  the  following  relation  for  determining 


V 


.'w  T{'“-  [“<^(¥3- -©■]}'“ 


S.  ,>,{-[•  -  ©1  ["&  -  '*]  -  ?} 

pa  f-z 


(232) 


l-l^] 

V  -1 


Va 


Vz. 

4 

*  t  ,  .  vl, 


2T 


paK  n  (p-v)  bt 


■!  Jf- 


KV(0-*O*bJ  J 


dz  *=  0. 


Equations  279  and  281  now  become 


?  sin^  *  Kp(/3-y)b1  , 


(293) 


f}*  cost  «  ±  \f  S4-K*fC31  (0-i/)*b*  , 


(294) 


which  substituted  in  Eq.  292  give 


& 


Ml^-Osin* 


f“  ■  ££}  [‘  *  (t)'  GS)]  r 


2Zj1paK 


If (r  +a)z  r 

— s-y 


Kr,  a)’ 
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32(ra+a)  z. 
5 (r  -a)4 


2a it 


+ 

fz‘ 

P*  sin<t> 

4a[l  -  f^Vl 

LKp (p-v) 

r  L  \*tJ  J 

fi*] 

_ 

2T 

'  *?- 

[•  +  pa/3* 

cos^  pa  1  - 

(295) 


x  ~  P*  ~  -  paXft(/3-p)tan^ldz  *  0. 

L  Va  J  z*  J 

Here  <p  replaces  bt  as  the  unknown  quantity,  the  two  being 
related  by  Eq.  293.  The  integral  in  Eq.  295  is 

P*sln0  »  »  M  -  f 2a(^->)*  _  „»]* 

W(/3-k)  p  *  15  zl  [  r,'a  ?  J 


+  palzi 


r  -a 

a 


iVsinft  f  2T 

V"  * 


-  paKp(£-»Otan<^  +  ("  ^a)  + 


+  2z 


]•’ 


-  paKp(0-iOtan$ 


(296) 


which  in  Eq.  295  gives 

P  sin^  +  Q  cos$  +  R  sin0  tan$  *  0 


(297) 


where 


P 


(P-v) 


4(r#+a) 

(r,-*)* 
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16(ra+a)z  r  4(ra+a) 

+  - —  ^paK  -  - — 


3(ra-a) 


1  + 


a  N2  I**"*' 


(ra-a)  L 


V  \v2¥a‘ 


+ 


32(ra+a)  zp  a 


+  L  P  a  16 

■  /  «  \  —  —  4LJ 


5(ra-a)4  J  15  1  L  ra“a 


lajl-JLl  _  £a 


+  pa  3  Zi 


ra”a 


+  2z 


1 


2T 


P  (  4? 


a  a  ,.a. 


"  2zlP  a-rKM(^) 


P  = 


2zx/3  n 
K 


(0-*O 


{}«{} 


3(ra-a) 


16(ra+a)‘ 

9(ra-a)4 


+  (r»+a)  As  .  16\ 

r.  -*•,*  Vs  9/ 


(ra-a) 


lo/3az 


15KP  (0 


Srf 


2  2 
a 


r  na 


-  5 


>pa[ 


+  2 ST*  +  Tl  f 15  25 

4z4  z4  V2  '  4 


-  2z  ^p2  a*  p2Kp  (p~v) 


2z3p 


P  = 


l'  ■  4(ra+a) 


(P-v)  if  paK  - 


4  *  (0  ©9] )' 


64 (r +a) 


45(ra-a) 


(rf-a) 

a+a)"l  ,  frf 2a(fi-r)*  ,al 

.a)4 J  15K M(P-v)  tVpa[  rt-a  '  fi  J 


5T 

2z‘ 


*T  +  -  2z  ^p2 a* p2Kp(P-v) 


16z yp  p 


>  Ar»+a)  f2  . 

b  +  ty  (v^JJ 
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x  I2  <V!L  [2  +  (V)1 

^  (ra-a)  L 


V* 


-  paK 


8(r  +a) 


45(ra-a) 


i6z  a 


15JQz(j8 


ra-a 


J  29!  /  47.  J 


(298) 


Q  =  T  P*aVz  -  tf2]  -  4pa02  (299) 

i  l  |_  ra-a  J  zx 


R  =  - 


p2a2/32z  I  — .  -  /321  +  4pa/S2  .  (300) 

1  L  I2  a  J  zi 


We  now  see  that 


R  -  -  Q  5 


(301) 


hence  Eq.  297  can  be  written 


P  sin0  +  Q(cos0-sin0tan0)  =  0  , 


—  —  12 

Psin0cos0+Q( cos  0-sin  < 


=  0  , 


P  sin20  +  Q  cos20  j  sec0  =  0  , 


(302) 


tan20 


-S. 


(303) 


Since  the  tangent  has  a  period  *  we  see  that  this  equa¬ 
tion  has  a  root  such  that 


0  <  20  <  t  , 
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or 


0  £  0  <  |  .  (304) 

Other  roots  of  Eq.  303  can  be  obtained  by  adding  multiples 

7 r 

of  7T  to  20 ,  or  multiples  of  to  0.  We  thus  obtain  four 

values  of  0  which  satisfy  Eq.  303  and  are  equally  spaced 
over  the  complete  circle  2tt.  Of  these  we  retain  only  those 
two  which  are  smaller  than  tt,  as  is  evident  from  the  figure 
on  page  12  and  the  fact  that  yQ  in  Eq.  39  must  be  negative. 

These  two  values  differ  by  and  of  them  we  wish  to  choose 

that  one  which  minimizes  Eq.  283,  or,  what  is  the  same  thing 
since  ~K  «  TJ,  T,  choose  that  for  which 


I 


+  c2 


)dz  =  minimum. 


(305) 


That  A  «  B,  C  follows  from  Eqs.  285,  286,  287,  and  291,  _ 

f rom_which  we  see  that  IT  is  of  second  order  in  £  >  whereas  B 
and  C  are  of  first  order.  Noting  Eq.  291  and  the  fact  that 
£  is  arbitrarily  small  but  fixed,  we  see  that  the  variable 
which  we  determine  in  minimizing  Eq.  305  is  bx ,  which  in 
view  of  Eq.  293  is  a  measure  of  0,  The  left-hand  side  of 
Eq.  297  and  hence  Eq.  302  was  obtained  by  dividing 


^  by  £ ,  thus 


1  dl  =  1  _dl 
£  db  dbx 


|  Psin20+Qcos20jsec0. 


(306) 


Noting  Eq.  293  it  follows  that 


_1_  d^I 
4*  db* 


_d_ 

d  0 


(  dl 

Vdbt 


dl\ 

J 


db 


l 

d0 


_  hiiLjlzJi )  -j  (Pcos20-2Qsin20)sec0 
P  COS0  ^ 


(307) 


+  ^  Psin20  +  §cos20jsec0tan0j  . 
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The  two  values  of  which  we  have  retained  satisfy  Eq.  302. 
For  each  of  these  Eq.  307  becomes 

2  2 

d_J_  =  (Pcos2$-2Qsin20) .  (308) 

db*  p  cos  <f> 

Since  p  >  v  and  K  is  positive  due  to  Eq.279,  the  sign  of 
this  expression  is  that  of 

P  cos20  -  2Q  sin2<£  .  (309) 

Since  our  two  values  of  0  differ  by  we  see  that  Eq.  309 

has  the  same  absolute  value  for  both,  but  opposite  signs. 

We  choose  that  value  of  <P  for  which  this  sign  is  positive, 
corresponding  to  a  minimum  of  I  (for  fixed  £),  and  discard 
the  other  value  of  <jf>,  which  corresponds  to  a  maximum. 

We  see  from  Eq.  303  that  the  two  possible  quadrants  of 
2 <(>,  and  hence  the  two  corrjes ponding  45-deg  sectors  of  0,  are 

determined  by  the  sign  of  For  each  of  these  the  signs  of 

sin2$  and  cos 20  are  also  determined.  It  follows  that  for 
each  of  these  the  sign  of  Eq.  3QJ3  is  determined  by  the  signs 
of  P  and  Q.  Given  the  signs  of  P  and  Q  we  can  hence  see  im¬ 
mediately  which  45-deg  sector  must  contain  in  order  that 
Eq.  303  be  satisfied  and  Eq.  309  be  positive.  These  sectors 
are  as  follows. 


Sign  of  P 

Sign  of  Q 

Sector  containing  ♦ 

4- 

- 

hk 

V 

V 

o 

- 

- 

A 

A 

toH 

- 

4 

1  <  0  < 

2  9  X 

4- 

4- 

4r  <  *  <  * 

In  any  particular  case  once  the  signs  of  ¥  and  have  been 
determined,  this  table  specifies  the  sector  which  contains 
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DETERMINATION  OF  THE  DAMPING  CONSTANT 
AND  PRE CESS ION AL  ANGULAR  VELOCITY 


Now  that_we  have  determined  <p,  b  follows  from  Eq.  279, 
after  which  x  follows  from  Eq.  93.  Noting  that  the  xy-axes 
of  the  figure  on  page  31,  have  their  origin  at  S  ana  are 
rotated  clockwise  by  an  angle  <f>  from  those  in  the  figure  on 
page  12,  whose  origin  is  at  W,  as  indicated  in  the  figure  on 
page  66,  we  see  that  x  and  y  on  page  12  are  given  by 


xrt  ~  4  +  x  co s<t>  -  I  + 


2ba  cos0 
3(r|-aa) 


^  +  2aaff*sin0cos0  V 

'  3(ra-aa)ty(/?-*V  * 


y0  = 


-  x  sin0  =  - 


2ba  sin <f> 

7/  a  a . 
3(ra-a  ) 


2a 

3 


aaflasina0 _ V 

(ra-aa)K*x(0-»O/ 


(310) 


These  are  the  coordinates  of  the  center  of  gravity  of  the 
mercury.  The  components  of  the  force  exerted  by  a  damper 
on  the  main  spinning  body  now  follow  from  Eq.  30,  thus 


*  ■  My0^a 

*  Mxo0*  = 


=  M(JJ  fi^jLkiiL± - X 

X3Cr*  -a'mUl-ry 
(i  + 

x  3(ra-aa  )Km  (P-vy 


(311) 


also  the  torque  components  exerted  by  a  pair  of  dampers  on 
the  main  spinning  body  are  given  by  Eqs .  3i,  34,  and  311, 
thus 


x  component  of  torque 
v  component  of  torque 
z  component  of  torque 


kx{- 

k  {, 


(312) 
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where 


kx  =  M/3  Lk  sin <f> , 


k  *  M/?  L(2  +  k  cos*)  , 


(313) 


k  =  M/3  k  sin <f>  =- 


4a3/3a  sin* _ 

3(ra-a2 )Km (P-v) 


(314) 


We  note  thav  the  coordinates  in  Eq.  310,  che  force  components 
in  Eq.  311,  and  the  x  and  y  torque  '•cmponen*  s  in  Eq.  312  are 
all,  closely,  proportional  to  4,  sii  e  in  _..ese  expressions 
the  coefficients  of  |  depend  upon  p  and  .  ,  which  vary  but 
little.  The  z  component  of  torque,  however,  varies  as  |a. 

From  Eq.  312  we  new  have 


r.  - 

;  =  ta 

'  A  s 


(315) 


hence  good  through  first-or  er  terms  in  £  ,  v  is  constant. 

The  results  given  by  the  approximate  picture  presented 
in  the  last  part  cf  the  section  between  pages  7  and  14  now 
follow  from  Eqs .  35  and  37,  which  give 


k  L 

i  =  -JL.  t 
*  Av  *  » 


(316) 


Reduction  in  preces- 
sionai  angular  velo- 
city 


(317) 


Here  Eq.  31/  can  be  used  to  determine  p.  Starting  with  the 

A  v 

approximate  value  p  &  we  determin '  the  right-hand  side  of 

0  Ay 

Eq.  317,  which  is  the  amount  by  which  must  be  reduced  to 
give  the  next  approximation  to  P.  Repeating  this  process 
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using  the  corrected  value  of  £  we  obtain  the  next  approxima¬ 
tion  to  0-~anu  =;o  on. 

In  Eq.  316  the  coefficient  of  £  is  essentially  constant; 
hence  £  is,  closely,  a  dying  exponential,  the  damping  factor 
(coefficient  of  t  in  the  exponential)  being 


a  =  - 


v.318) 


The  negative  reciprocal  of  a  is  the  time  constant. 

The  results  given  by  the  approximate  picture  presented 
previously  (pp.  25-28)  follow  from  Eqs.  73  and  74,  thus, 
noting  Eq.  312, 

.  k  L 

dt  'tP)  =  £> 


or,  if  /3  is  considered  to  be  essentially  constant. 


*  “  *  3/J  ^  * 


(319) 


Also 


Y  (Ar-B/S)  -  k  £, 


Reduction  in  preces- 
sional  angular  velo¬ 
city 


h  k  L 

Av  0  v 

T  "  * 


(32(  ) 


The  dampir.g  facto"  corresponding  to  Eq.  319  is 

k  L 

“  ■  -  ~tp- 


(321) 


Equation  317  differs  from  Eq.  20,  and  Eq.  318  differs  from 
Eq.  321  in  that  Av  appears  in  place  of  B£;  however,  the  ratio 
of  these  two  quantities  shoula  approximate  unity,  since 
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-g-  approximates  0.  It  follows  that  the  results  given  by  the 

two  pictures  are  m  approximate  agreement.  It  may  be  noted 
that  in  Eqs.  319,  320,  and  321  0  is  the  angular  velocity 
seen  by  an  observer  who  is  located  at  the  instantaneous 
center  of  curvature  of  the  path  of  S. 


APPLICATION  OF  PROCEDURES  FOR  LINEAR 
DIFFERENTIAL  EQUATIONS 


Since  the  coefficients  of  |  in  the  first  two  equations 
of  Eq.  312  are  closely  constant  it  appears  that  we  can  check 
and  perhaps  improve  the  treatment  in  the  section  on  pages  25- 
28  by  using  the  methods  available  for  solving  systems  of 
linear  differential  equations.  Let  us  choose  a  stationary 
set  of  XY  axes  in  the  plane  of  the  upper  damper  with  its 
origin  at  W,  in  the  vicinity  of  S,  as  shown  in  the  figure 
below.  £  -  WS,  and  is  hence  the  x  coordinate  of  S.  Knowing 


the  torque  components  aJ one  the  x  and  y  axes  from  Eq.  312, 
we  can  obtain  those  along  the  X  and  Y  axes  immediately,  thus 


X  component  of  torque  =  cos 0  -  k^4  sin0 

-  v  -  v- 

Y  component  of  torque  *  k  t  sinO  +  kj  cos0 

x  y 

*  k  Y  +  k  X  . 
x  y 


(322) 
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The  differential  equations  for  the  motion  of  S  are  hence, 
closely, 


R  •»  fly  • 

—  x  +  “  v  ~  k  Y  +  k  X 
L  L  *  x  y 

fly  ’  B  •• 

“  r  x  —  y  s  y  x  -  V:  Y 
I.  L  x  "  y  * 


(323) 


or,  denoting  by  D, 


(BD*  -  Lky)X  +  (Aud  -  Lkx)Y  -  0 

*  a  (324) 

(AvL  -  Lk  )X  -  (BD-  Lk  )Y  =  0. 

x  y 


Equating  the  determinant  of  the  system  to  zero  we  obtain 

(BD*-Lk  )  (AyD-Lk  ) 

y  ,  x  -  0,  (325) 

(AyD-Lk  )  -(BD-Lk  ) 

x  y 


-(3D*-Lk  )*  -  (AyD«Lk  )*  -  0, 

y  ^ 

BD*  -  Lk  *  *  i( AyD-Lk  ).  (326) 

y  * 

Equation  325  is  a  quartic  polynomial  in  D,  whose  four  roots 
are  composed  of  the  two  pairs  of  roots  of  Eq,  326  that  are 
obtained  for  the  two  possible  signs.  Since  it  is  evident 
that  the  conjugates  of  the  two  roots  obtained  for  one  sign 
in  Eq .  326  are  roots  of  this  equation  with  the  opposite 
choice  of  sign,  it  follows  that  we  need  consider  but  one 
sign.  Choosing  the  plus  sign  Eq,  326  may  be  written 

BD*  -  iAyD  -  L(k  -ik  )  *  0;  (327) 

y  ^ 


hence 
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iAi>  Ti  .  /  "i  43L  /i.  ,  \ 

— r^r  1  i  s/  1  -  — r — r  Ik  ) 

2B  J_  id  y  x' 


(328) 


For  either  choice  of  sign  we  obtain  a  root,  whose  conjugate 
is  also  a  root  of  Eq.  32S,  thus 


D  =  a  Jr.  i/3. 


(329) 


Corresponding  to  such  a  pair  of  roots  we  have  a  mode,  for 
which  che  damping  constant  is  a,  and  the  precessional  angular 
velocity  is  0(0  >0). 

a  and  0  can  be  obtained  exactly  from  Eq.  328;  however, 
by  applying  the  binomial  theorem  to  the  radical  in  Eq.  328 
we  can  obtain  the  following  approximate  expression 

i  _i 

ri  iAi/  f.,  .  r,a  ,  1  ,’i  /  4BL  „ 

D  *  ^1±|1  +  -r  1  -  ~  (k  -  lk  )  ) 

2B  i  t  2  \  .2  a  y  x  J 


.  l  (.  l\  i_i  /  _1BL  ,,  ,Y  ]  \ 

+  2  {  2J  2  V'  AV  (ky"lkx V  +  ‘"J  J 
%  h  ±  (k  ,ik  , 

23  i  t  a*j  y  x 


- -  (k  -k  -i2k  k  )  -t  .  . . 

.4  4  v  x  x  y 
A  v  1 


(330) 


Noting  Eq,  329  this  gives  for  the  plus  sign 


Lk  23L  k  kw 
r,  X  X  y  , 

D  “  -  -XT  -  ~TT  +  '  ■  ■ 

A  v 

4  A*i  _  — y  .  (v*-k*)  + 

LB  »*„»  y  * 


(331) 


whence 


o  *  - 


2BL  k  k 


.j  s 

A  v 


+  .  . 


(332) 
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R  =  A t 

p  B 


Lky 

Av 


BL 


*  3 
A  v 


-  (k*-k3) 
3  v  y  X 


We  note  that  through  first-order  terms  in  k  and  k()  Eq.  3.32 
agrees  with  Eqs.  317  and  319.  ” 


The  corresponding  mode  can  be  obtained  us_ing  complex 
notation  for  damped  sinusoids;  thus  if  X  and  Y  denote  the 
complex  numbers  corresponding  to  X  and  Y,  respectively,  and 
D  =  a  +  ifi  the  differential  equations  (Eq.  324)  become* 


( BD* -Lk  )X  4  ( ArD-Lk  )Y  =  0 
y  x 

(AI/D-Lkx)X  -  (BD2-Lky)Y  =  0. 


(333) 


The  solution  of  this  set  or  linear,  homogeneous,  algebraic 
equations  can  be  written  down  immediately  by  noting  that  X 
and  Y  are  proportional  to  the  cofactors  of  the  elements  of 
either  row  of  the  determinant  in  Eq.  325--say  the  first, 
thus 


X  -  -  e  (BD* -Lk  ) , 

y  (334) 

Y  -  -  c  ( Ai'D-Lk  ) 

A 


7  Complex  notation  for  damp  d  sinusoids  can  easily  be  ob¬ 
tained  as  foil  owe.  Let  a  one-to-one  correspondence  between 
damped  sinusoids  of  fixed  damping  constant  a  and  angular  fre¬ 
quency  cj,  and  points  of  the  complex  plane  be  defined  by  the 
relation  ^ 

Ac  sin(wt40)  ~  A  /_  $  *  A(cos<H-isinG) ; 

then  it  can  readily  be  shown  that  the  complex  number  corre¬ 
sponding  to  a  linear  combination  of  damped  sinusoids,  all 
having  the  same  a  and  is  equal  to  that  same  linear  com¬ 
bination  of  the  complex  numbers  corresponding  to  the  compo¬ 
nent  damped  sinusoids;  also  that  the  complex  number  corre¬ 
sponding  to  the  time  derivative  of  a  damped  sinusoid  is 
equal  to  (cH-ikO  times  the  complex  number  corresponding  to 
the  damped  sinusoid  itself. 
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where  cx  is  an  arbitrary  complex  constant.  Since  Eq.  326 
with  the  plus  sign  is  valid  Eq.  33<"  becomes 


X  =  ic, 
Y  =  c 


(335) 


where  c  =  |c[  /  0  is  an  arbitrary  complex  constant.  It 
follows  that 


X  =  |  c  |  e0*  cos(/3t+0c) , 
Y  =  1  c  |  e0^  sin(/3t+0c); 


(336) 


also 


|  -  yp+Y*  =  | c  |  e0*. 


(33?) 


The  phase  of  the  damped  sinusoids  Eq.  336  is  hence  such  that 
S  spirals  around  W  with  an  exponentially  decreasing  radius. 


It  might  appear  that  another  mode  can  be  obtained  by 
using  the  minus  sign  in  Eq.  328,  and  that  the  motion  of  S  is 
a  composite  of  the  two  modes;  however,  for  the  torque  equa¬ 
tions  (Eq.  312)  to  be  valid  it  is  necessary  that  the  origin 
be  at  W,  on  the  precession  axis.  Furthermore  with  the  minus 

Lk 


sign  the  value  of  p  obtained  approximates  —j^j,  as  can  be  seen 

from  Eq.  330.  At  this  extremely  low  precessional  angular 
velocity  the  above  analysis  of  the  action  of  the  mercury 
damper,  upon  which  Eq.  512  is  based,  is  no  longer  valid. 


POLAR  FORM  OF  THE  TORQUE  DIFFERENTIAL  EQUATIONS 

Referring  to  the  figure  on  page  98  we  note  in  regard  to 
the  motion  of  S  tnat 

Radial  component  of  velocity  ■  4, 

Transverse  component  of  velocity  *  $0.  (338 

M  |  |  ' 

Radial  component  of  acceleration  *  4-4®  * 

•  •  •• 

Transverse  component  of  acceleration  »  240+40. 
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Noting  that  p  -  Q  it  follows  that 

X  «  -  L  (2«+*«  “  2tx 

E  +X  to  -  2V 

closely  where  2 tx  and  2ty  are  the  x  and  y  components,  respec¬ 
tively,  of  the  total  torque  exerted  by  both  dampers  on  the 
main  spinning  body.  These  relations  may  be  written 


(Au-2Bp)i  -  B| p  -  2Lrx 
B  '£  +  (Av-Bp)pt  =  2Lry. 


(339) 


Applying  Eq.  339  to  the  case  just  considered  in  the 
previous  section  we  have,  noting  Eqs.  312  and  337, 

i  -  t-a'i,  0-0,  2rx=kx4,  2ry  =  ky{;  (340) 


hence 


( Av-2Bp)oi  -  Lk 

x 

Ba*  +  (Av-Bp)p  ~  Lky. 


These  relations  may  be  written 


Av 

B 


Lk 


a  = 


Av- 2Bpy 


-  P 


Lk^ 
B  P 


(341) 


(  342) 
(343) 


which  give  the  damping  constant  and  the  reduction  in  the  pre- 

Am 

cessional  angular  velocity.  Again  we  note  that  since  -g- 

approximates  0,  Eq.  342  is  in  approximate  agreement  with  Eqs. 
318  and  321.  In  using  Eqs.  342  and  343  p  is  first  assigned 

Av 

the  approximate  value  -g-,  and  a  first  approximation  to  Of 
calculated  from  Eq.  342.  This  and  the  approximate  value  of 
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/3  placed  in  the  right-hand  side  of  Eq.  343  give  the  amount 
Ay 

by  which  -g~  must  be  reduced  in  order  to  obtain  the  second 

•  A  v 

approximation  to  p.  Using  this  instead  of  -g-  we  can  now 

repeat  the  above  process  to  get  a  second  approximation  to  a, 
a  third  approximation  to  0,  and  so  on. 

Finally  we  note  that  the  equations  (Eq.  339)  provide  a 
better  means  of  determining  damping  behavior  and  precessional 
angular  velocity  than  does  the  method  described  in  the  sec¬ 
tion  on  pages  25-28.  We  shall  see  that  these  equations  will 
prove  useful  in  considering  the  case  of  the  crescent  configu¬ 
ration  of  the  mercury,  for  which  the  torque  components  are 
given  by  expressions  which  differ  greauly  from  Eq.  312,  and 
for  which  therefore  linear  differential  equation  procedures 
cannot  be  used. 


CASE  WHERE  THE  MERCURY  CONFIGURATION 
IS  CRESCENT- SHAPED 


We  shall  next  consider  the  case  where  the  mercury  is 
crescent- shaped,  and  completely  tills  an  angle  2flm  of  the 
channel,  or  race,  the  flow  being  turbulent.  In  so  doing  we 
shall  use  che  same  notation  which  we  have  been  using  in  con¬ 
nection  with  the  eccentric  angular  configuration.  Instead 
of  the  figure  on  page  31  we  now  have  the  configuration  shown 
in  the  figure  on  the  following  page.  The  volume  of  the  mer¬ 
cury  is  evidently 


Volume  of 
mercury 


2zi*mCVri)*  (544) 


Again  denoting  the  x  coordinate  of  the  center  of  gravity  of 
the  mercur>  by  x  we  have 


x  times  mercury 
volume 


r*co5$d0drdz 
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Here  0m  is  determined  by  Eq.  344,  since  the  mercury  volume 
is  known. 

If,  as  before,  we  observe  the  system  while  riding  with 
axes  which  rotate  counterclockwise  at  the  precessional  angular 
velocity  /?  we  see  the  configuration  shown  in  the  above  figure. 
Noting  the  figure  on  page  12  we  now  take  moments  about  S,  as 
described  in  the  section  on  pages  35  and  16,  and  thus  obtain 

lFi  U  =  047) 

Here  rm  is  the  mean  radjus  of  the  mercury,  and  is  the 

velocity  of  the  mercury  relative  to  the  channel.  Kq  is  a  co~ 
eificient  obtained  from  Chezy’s  formula,  which  coefficient 
depends  upon  the  shape  of  the  channel,  the  area  of  the  surface 
of  contact  of  the  mercury  with  the  channel,  and  the  density 
and  viscosity  of  the  mercury.8 


Rouse,  Hunter,  and  J.  W.  Howe.  Basic  Mechanics  of 
Fluids.  New  York,  Wiley,  1953.  P.  149. 
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Noting  Eqs.  32,  34,  and  the  fact  that  there  are  two 
dampers  we  obtain  the  following  components  of  the  torque 
exerted  by  the  mercury  on  the  main  spinning  body. 


Total  x  component  of  torque  = 


2r  =  2 |F  |L 
x  1 

2LK  r3(/3-v)2 
o  m 


(348) 


Total  z  component  of  torque  -2 t  =  2  j F  |£ 

Z  X 


*  2kV(^-02. 

o  m 


(349) 


Equation  347  in  the  first  equation  of  Eq.  30  now  gives 


y0  =  - 


K  rS (/3-v)2 

o  _ _ 


(350) 


whence,  noting  Eq.  346  and  the  figure  on  page  106, 


sin <f>  =  -  ^r> 


(351) 


which  relation  determines  <p. 
106  we  now  have 


Again  noting  the  figure  on  page 


x0  «  $  +  x  COS0, 


(352) 


whence  the  second  relation  in  Eq.  30  gives 


IF, I  -  M01 (4+xcos0) , 


(353) 


It  now  follows  from  the  second  equation  of  Eq.  34  that  for 
the  two  dampers  together 

Total  y  component  of  torque  «  2r„  ■*  2|F„|L 

y  *  (554) 

-  2UW*(t-hceos^). 


•mv 

•  W  Hill  1^  ! 

pf 
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Equating  Eq.  34°  co  A v  we  obtain  for  the  rate  of  increase 
of  the  angular  velocity  of  spin 


v 


2K„rh/J-K)a 
o  m _ 

A 


(355) 


Finally  let  us  apply  Eq.  339  to  the  present  case  in 
order  to  obtain  information  on  the  damping  and  the  decrease 
in  precessional  angular  velocity.  In  so  doing  we  shall 
neglect  the  variation  of  p  and  v.  Substityting  Eq.  348  in 
the  first  equation  of  Eq.  339  and  placing  £-0  we  obtain 


.  2L*K  r*(/3-n)2 

(Aj/-2B/?)£  =  - - .  (  356) 


Letting  £  denote  the  initial 
and  integrate,  thus 


(Av-2Bp) 


value  of  £  we  separate  variables 


-  2l’k 

o  m 


or 


I 


4L2K r*(0-r)a 
_ o  m _ 

2B/3-  Av 


t . 


(357) 


Finally  we  substitute  Eq.  354  in  the  second  equation  of 

~B  “  Pj' 

The  decrease  in  the  precessional  angular  velocity  is  thus 
found  to  be 

T  ~  P  m  ifa  [2lWU  tIcos0)-B*].  (358) 


•  • 

Hero  it  is  likely  that  in  the  bracket  £  and  £  can  be  neglected 
but  this  can  best  be  decided  in  specific  numerical  cases.  In 
this  connection  we  no’re  that 
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*  =  “  Js  (  2B8^Aif  J  ’  (359) 

which  can  be  obtained  by  differentiating  Eq.  357  twice. 

Several  things  are  immediately  evident  from  the  above 
analysis.  First,  we  see  from  Eqs.  '50  and  351  that  as  4  de¬ 
creases  during  the  course  of  the  damping  (-yp)  and  £  increase. 
The  smallest  value  of  4  that  is  compatible  with  the  crescent 

configuration  is  that  for  which  -y  =  x,  in  which  case  <t>  = 

For  smaller  values  of  4  the  mercury  configuration  becomes  ec¬ 
centric  annular. 

Second,  we  note  from  Eq.  357  that  4  decreases  parabolli- 
cally  instead  of  exponentially,  as  shown  below.  The  rate  of 


damping  increases  as  4  decreases,  and  is  greatest  just  before 
the  crest ent  flow  breaks  and  goes  over  into  the  eccentric 
annular.  In  fact  differentiating  Eq.  .357  we  obtain 


% 

4 


1 

4  \  2B0-A* 


(360) 
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from  which  we  see  that  the  rate  of  damping  varies  inversely 
as  closely,  instead  of  directly  as  which  is  the  case 
with  the  eccentric  annular  configuration. 

Finally,  we  see  that  since  the  first  part  of  the  curve 
of  4  against  t  is  a  parabola  which  is  concave  downward,-  and 
the  last  part  of  this  curve  is  a  damped  exponential,  which 
is  concave  upward;  and  since  there  is  an  intermediate  part 
of  this  curve  which  corresponds  to  values  of  £  for  which  the 
flow  is  eccentric  annular  but  very  eccentric;  it  follows 
that  the  £t  curve  as  a  whole  has,  if  differentiable,  a  point 
of  inflection,  and  in  any  case  should  be  capable  of  being  ap 
prcximated  by  a  straight  line  over  an  appreciable  range. 

Such  a  straight  line  is  shown  dotted  in  the  figure  on  page 
109. 
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Appendix  A 

EFFECT  OF  THE  TILT  AND  THREE  DIMENSIONAL  NATURE 
OF  THE  DAMPER  ON  THE  FORCE  WHICH  IT  EXERTS 
ON  THE  MAIN  SPINNING  BODY 


The  force  exerted  by  the  mercury  on  the  main  spinning 
body  was  derived  on  pages  8-12.  In  this  derivation  the 
damper  was  considered  to  be  essentially  two  dimensional,  its 
plane  being  fixed,  and  its  motion  being  in  this  plane.  This 
simplified  the  derivation,  and  gave  for  the  resultant  force 
exerted  by  the  mercury  on  the  main  spinning  body  a  force 
which  lies  in  the  plane  of  the  damper.  Actually  the  damper 
is  three  dimensional,  and  is  tilted  slightly.  This  gives 
rise  to  a  torque,  which  may  shift  the  resultant  force  away 
from  W,  and  which  may,  in  addition,  give  rise  to  a  torque 
whose  vector  is  parallel  to  that  of  the  resultant  force. 

The  purpose  of  this  appendix  is  to  refine  the  above  deriva¬ 
tion,  taking  into  account  the  tilt  and  three  dimensional 
nature  of  the  damper,  and  see  whether  or  not  the  results  ob¬ 
tained  differ  appreciably  from  those  obtained  above. 

Determination  of  FT .  Let  us  place  a  set  of  xyz  axes  with 

the  z  axis  coincident  with  the  precession  axis,  the  +z  direc¬ 
tion  being  outward  away  from  the  body,  and  the  positive  x 
axis  passing  through  the  center  of  gravity  of  the  mercury. 
Noting  Eq.  6  it  follows  that 


0l  (ix+jy)dm  =  /3a 


(*t 


x  dm  +  j 


Fl  =  i/3aMx 


(»1) 


where  the  integrations  extend  throughout  the  mass  of  the 
mercury,  and  x  is  the  x  coordinate  of  the  center  of  gravity 
of  the  mercury.  We  thus  see  that  the  total  centrifugal  force 
is  identical  with  that  obtained  by  concentrating  all  mass  at 
the  center  of  gravity.  This  re cult  is  in  agreement  with  Eq. 
14. 
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Determination  of  F^.  Equations  15  and  16  are  formally 

valid  in  the  present  case;  hence  we  again  obtain  the  result 
(Eq.  17),  namely 


Fc  -  0.  ( A2) 

Since  the  total  Coriolis  force  is  zero  the  total  force  -F  is 
composed  of  F^  alone,  as  before. 

Determination  of  r^.  Let  us  place  a  set  of  xyz  axes  with 

the  z  axis  coincident  with  the  precession  axis,  the  positive 
z  direction  being  outward,  and  the  positive  x  axis  passing 
through  S,  as  shown  in  figure  below. 


The  torque  vector  due  to  the  centrifugal  force  field,  calcu¬ 
lated  for  the  origin  0  in  then,  noting  Eq.  6, 


T 


L 


/32  (ix+jy+kz)  x  (ix+jy)dm 


i  j  k 

xyz 
x  y  0 


dm 


v.  u* 
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i  \  yzdm  +  j 
JM 

Next  let  us  place  a  set  of  x,y1z1  axes  with  the  zx  axis 
coincident  with  the  spin  axis,  the  positive  direction  being 
outward,  the  origin  at  S  and  the  positive  xt  axis  passing 
through  the  center  of  gravity  of  the  mercury  (see  figure  on 
previous  page  and  the  figure  below).  The  xyz  and  xyz 


xzdrn^). 


(A3) 


coordinates  of  any  point  are  then  related  by  the  equations 
x  -  xl  cos<p  cosa  +  yx  sin*  cosa  +  z1  sinat  +  L  sina 
y  -  -  x^  sin  <p  +  y1cos0  (A4) 

z  =  -  xJL  cos<£  sina  -  y^  s±n<P  sina  +  z ^  cosa. 

Substituting  Eq.  A4  in  Eq.  A3  we  obtain 

=  -  i£*  j  [x*  sin</>  cos0  sina 

M 

.  2  2 

-  y*  sin0  cos0  sina  -  x^  (cos  <t>-s±n  <p) sina 
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-  xJz1  sin #  cosa  +  y1zx  cos#  cosa]dm 

+  j  P*  \  f-  x*  cos a#  sina  cosa 

1 

2  2  2 

-  yi  sin  #  sina  cosa  +  z  sina  cosa 

-  2  xJy1  sin#  cos#  sina  cosa 

+  xjzl  cos#  (cosaa-sinaa)  +  y^  sin#  (cos*a 

2  2  2 

-  sin  a)  -  xiL  cos#  sin  a  -  yiL  sin#  sin  a 

+  sina  cos  a]  dm. 


Since  the  x.yj  plane  and  the  XjZL  plane  are  planes  of  sym¬ 
metry  for  the  mercury  the  various  products  of  inertia  vanish, 
leaving 


=  -  i  p2  sin#  cos#  sina  \  (x*-y*)dm 

JM  1 

+  j  0*  sin2a  E-(xJ+yJ)  (yJ-x*)cos2# 
+  2za ]dm  -  MLxJ  cos#  sin*aj- 


(A5) 


where  xx  is  th£  xt  coordinate  of  the  center  of  gravity  of 
the  mercury.  yx  and  zx  are  zero  because  of  symmetry.  We 
shall  next  digress  long  enough  to  compute  the  relevant 
integrals  using  the  procedure  and  cylindrical  coordinate 
system  already  used  in  Eqs .  87  and  90.  In  so  doing  we  shall 
omit  terms  of  order  higher  than  one  in  b. 


We  have 


)dm  *  p 


ra (cosa  6 

COS0 


-  sin* 0)rdzd0dr 
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“  *  LX  f  'J- Wb[l  -  (f- 


an 


=  o, 


cos0jcos20dzd0 


(A6) 


z  r 

i  n  nr*  a 


+-  ©‘ 


r  dzdfldr 


cose 


z 

a  4  L  {ri-«4+»a’b[l  -  (f-')  ]cos6}dzda 


RJL 

2 


X  (r»4-a‘ 


)dz 


=  P»  (rj-a4), 


z  r 

l  a*/*  a 


L 1 1 


z  rdzdfldr 


cos0 


a  2  L  L  2,{r5’a‘+2ab[i"  X)  ]c°s<’}dzd'’ 


dz 


z 

,  a  f  1  * 

=  pir(r  -a  )  \  z 

2  /I  2  v  3 

“  J  "»<Va  >V 

Equations  A6,  A7,  and  A8  in  A5  now  give 

tl  a  j  p*  ^2  sin2a  pirCrJ-a1  )zt  (4zJ-3rj~3a*) 

-  MLXj  cos0  sin*aj. 


(A7) 


(A8) 


(A9) 
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From  Eq.  291  we  see  that  b  ana  £  are  of  the  same  order; 
hence  since 


£  =  L  sina  (A10) 

it  follows  that  |,  b,  and  sina  are  of  the  same  order,  and 
the  last  term  in  the  bracket  of  Eq.  A9  can  be  omitted.  Sub¬ 
stituting  from  Eq.  89  we  then  obtain  finally 


r,  a  j  sin2a  (4z?-3r*-3a2).  (All) 

Li  *4  12 

Determination  of  Tq.  First  we  note  that  since  Fc=  0,  as 

indicated  by  Eq.  A2,  the  torque  vector  due  to  the  Coriolis 
force  field  is  independent  of  the  point  for  which  it  is  cal¬ 
culated.  We  may  therefore  choose  S  as  the  point  for  which 
we  calculate  tq.  Next  we  note  that  since  the  cross  product 
is  distributive  tq  may  be  regarded  as  the  sum  of  the  contri¬ 
butions  of  the  f ollowing_£omponents  into  which  the  preces- 
sional  angular  velocity  p  is  resolved. 


1. 

P 

cos  a 

in  the 

+z 

l 

direction, 

2. 

P 

sina 

cos 4>  in 

the 

-X 

l 

direction , 

(A12) 

3. 

P 

sina 

sin0  in 

the 

■y. 

direction . 

We  shall  now  compute  these  contributions. 

The  Coriolis  forces  on  two  symmetrically  placed  equal 
elements  of  mass  dm  due  to  the  first  of  the  components  (Eq. 
A12)  are  shown  in  the  figure  at  the  top  of  the  following  page. 
It  is  now  evident  that  the  four  Coriolis  forces  which  act  on 
four  equal  elements  of  mass  dm  symmetrically  located  at  points 
(x.,±yi )  together  have  a  resultant  which  passes  through  S, 
and  hence  contributes  nothing  to  the  torque  vector  for  this 
point.  Since  the  mercury  may  be  regarded  as  composed  of  such 
sets  of  symmetrically  placed  mass  elements  it  follows  that 
the  contribution  of  the  first  (and  largest)  of  the  three 
angular  velocity  components  (Eq.  A12)  to  is  zero. 

Turning  to  the  second  of  the  components  in  Eq.  A12  we 
have  the  situation  shown  in  the  figure  at  the  bottom  of  the 
following  page,  where  Q  indicates  "up"  and  0  indicates  "down”. 
The  resultant  of  these  forces  is  a  couple  whose  vector  points 
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in  the  +y1  direction  and  is  the  desired  contribution  to  tq  . 

We  wish  to  determine  this  accurate  through  first-order  terms 
in  b.  But  the  angular  velocity  component  has  a  factor  sina, 
which,  we  saw  above,  is  of  first  order  in  b;  hence  ir.  comput¬ 
ing  the  moment  of  the  Coriolis  force  system  about  the  yx  axis 
we  may  place  b=0  in  the  expressions  which  give  and  the 
shape  of  the  mercury.  In  this  way  we  obtain  directly  the 
desired  first  term  in  the  Maclaurin  series  expansion  of  the 
moment.  Under  these  circumstances  we  see  the  mercury  as  an 
annular  ring  with  no  eccentricity  rotating  clockwise  with 
angular  velocity  ( >3- ^ ) »  as  shown  below. 


Using  cylindrical  coordinates  for  the  integration  we  obtain 
for  the  desired  moment 


z  r 

l  p  l  w  n  2 


f  ir  r 

Moment  »  2p/?sin«costf>  \  \  \  (^-OrcosflrcosOrdzdfldr 

*  Z  ^  0  cL 

*  4  p(3((3-v)  sinacos0(r^-a4 )*2z 
Moment  a  —  M/S( (i-v)  (rj+aa  )cos</>  sina. 


( Ai  3) 
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The  contribution  of  the  second  component  in  Eq.  A12  to  Tp  is 
a  vector  having  the  +y  direction  and  the  magnitude  (Eq,  A13). 

The  third  component  in  Eq.  A12  can  be  treated  in  the 
same  manner  as  the  second.  Since  the  two  are  perpendicular, 
and  one  has  a  factor  sin<£  whereas  the  ether  has  a  cos$>,  cor¬ 
responding  modifications  must  be  made  in  the  contribution  to 
ro  Noting  Eq.  A13  it  follows  that  the  contribution  of  the 
third  component  of  Eq.  A12  to  is  a  vector  in  the  -x  direc¬ 
tion  of  magnitude 


Magnitude  s  |  Mj3(i3-p)(r|+a,)sin^  sina.  (A14) 

If  these  two  vector  contributions  are  resolved  into  cowpo 
nents  along  the  x,  y,  and  z  directions  and  combined,  we  find 
that  tq  consists  of  a  vector  extending  in  the  4y  direction 
and  given  by 


j  M/3(/?-y)  (r*+aa  )sina. 


(A15) 


In  view  of  Eqs.  A1  and  A2  we  see  that  the  expressions 
in  Eq.  311  for  1 Fx |  and  1f*|  remain  valid  in  the  present  case. 
We  note  that  these  are  of  first  order  in  4  and  hence  b.  It 
is  for  this  reason  that  expressions  accurate  through  first 
order-terms  in  b  are  adequate  in  computing  torque  vectors. 

The  torque  vectors  tj,  and  tq,  given  by  Eqs.  All  and  A15, 
respectively,  are  evidently  provided  for  by  letting  the  force 
component  vector  Fl  act  at  W  as  before;  but  moving  the  point 
of  application  of  the  component  vector  Fj  to  a  point  on  the 
precession  axis  whose  distance  from  W  is* 


Distance  from  „  j  component  of  Cy^) 

Fa  to  w  a  ffpr 

i  >  "l 

a  _i_  A  +  2a* /r  sinftcosfl  \  (A16) 

2L  5(r|-a*)X*i(0-i'K 

x  {i  (3r*+J«*-«*J)  -  fl  -  ■§)  (r|+a*)} 

toward  the  center  of  gravity  of  the  main  spinning  body.  Here 
|F  !  is  given  by  Eq.  311. 


s 
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We  have  seen  that  the  torque  on  the  main  spinning  oody 
due  to  is  of  importance  in  determining  the  damping  factor 
or  and  the  rate  of  increase  of  the  spin  velocity  *,  whereas 
the  torque  due  to  Fa  is  of  importance  in  determining  the 
shift  in  the  precessional  angular  velocity  /3.  We  have  not 
altered  F . ,  and  have  changed  only  the  point  of  application 
of  Fj .  The  following  final  results  are  now  evident.  The 
three  dimensional  nature  of  the  damper  and  its  tilt  require 
no  correction  whatever  j n  the  damping  factor  a  or  the  rate 
of  increase  in  spin  velocity  T.  The  shift  in  prece s s ional 
angular  velocity  is  affected  by  this  refinement,  and  can  be 
gotten  by  replacing  L  by  [L-distance  given  by  Eq.  A.16]  in 
the  expression  for  ky  in  Eq.  313  before  using  it  in  Eq.  317. 

We  are  really  interested,  however,  only  in  the  damping  factor. 
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Appendix  B 

OUTLINE  OF  THE  PROCEDURE  FOR  SOLVING  A  PROBLEM 

A .  Eccentric  Annular  Configuration 

1.  Determine  a,  the  inner  radius  of  the  mercury 
with  no  wobble.  This  can  be  done  using  Eq. 
89. 


2 

T 

J  , 


4 

5 


Determine  K  using  Eq.  280. 

Tentatively  placing 

_  hiL 

~  IT 

determine  P  from  Eq.  298,  and  Q  from  Eq.  299. 

Determine  <f>  from  Eq.  303  and  the  table  on 

pagt  94. 

Determine  k  from  Eq.  314;  then  determine 

k  ,  k  ,  and  k  from  Eq.  313. 
x  y  z 


Obtain  v  from  Eq.  315 
and  a  from  Eq.  313. 


If  desired  the  calculated  value  of 


^  a  \ 

T  '  “) 

ive  an  improved  value  of  ,8,  which  can  then  be  used  to  replace 
q.  Bl;  and  the  subsequent  calculations  repeated,  and  so  on. 


(Bl) 


from  Eq.  317, 

can  be  used  to 


B .  Crescent  Configuration 


1. 

2, 


3. 


Determine  0m  from  Eq.  344. 

Determine  rm,  the  mean  radius,  from  the  relation 


m 


i  <v  v  • 


(  B2) 


Determine  X0  from  Chezy’s  formula,  which  is  de¬ 
scribed  in  the  reference  (footnote  8)  on  page 
106,  and  in  other  books  on  applied  hydrodynamics. 

Tentatively  placing 


T 


B3) 
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v  is  given  by  Eq.  355,  the  motion  of  the  spin 

axis  is  given  by  Eq.  357,  and  is  given 

by  Eq.  358,  in  which  £  and  |  are  given  by  Eqs. 
357  and  359,  respectively. 


a 


which  is  now  a 


Av 

If  desired  the  calculated  value  of  ( — g-  -  /9 

function  of  t,  can  be  used  to  suggest  a  value  of  (3  which 
supersedes  that  given  by  Eq.  B3,  after  which  the  subsequent 
calculations  are  repeated;  and  so  on. 


The  time  at  which  the  wobble  would  disappear  if  the 
crescent  configuration  did  not  go  over  into  the  eccentric 
annular  configuration  is  given  by  placing  |=0  in  Eq.  357  and 
solving  for  t.  Actually,  as  stated  on  page  109,  the  crescent 
configuration  cannot  persist  for  values  of  £  which  are  ^mailer 
than  that  given  by  equating  -  y0,  given  by  Eq.  350  and  x, 
given  by  Eq.  346.  The  value  of  t  at  which  this  value  of  £  is 
reached  is  given  by  Eq.  357. 


NOTE.  In  this  report  the  analysis  of  the  eccentric  an¬ 
nular  case  is  intended  for  use  only  in  connection  with  large 
mercury  dampers,  wherein  surface  tension  plays  but  a  small 
role.  (See  footnote  5,  page  79.)  In  any  particular  case  it 
is  hence  necessary  that  the  terms  in  the  calculations  which 
are  due  to  surface  tension  be  relatively  small. 
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